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Abstract 

We would like to discuss, within the inherent historical context, 
upon one of the main works of fundamental physics which has led to the 
formulation of the early resonance dual models prior to string theory. 

In fact, we shall focus on certain aspects of the fundamental Tullio 
Regge work of 1959, to be precise, on certain mathematical methods 
handled by him to pursue his original intentions mainly motivated 
to prove some double dispersion relations into the framework of non- 
relativistic potential scattering. 

It is well-known that the celebrated Gabriele Veneziano work of 1968 is al¬ 
most unanimously considered as the first step toward the formulation of 
resonance dual models of particle physics, from which then string theory 
soon arose. Nevertheless that, we are of the historiographical opinion that 
the early prolegomena to the first resonance dual models is the famous 1959 
seminal paper of Tnllio Regge, in which, in introducing complex orbital mo- 
mentua to determine further possible connections between potentials and 
scattering amplitudes, a notable formal technique was introduced, namely 
the so-called Watson-Sommerfeld transforrr^. Thenceforth, such a formal 
tool even more played a fundamental and pivotal role in theoretical particle 
physics, so that we would like to highlight the chief historical moments which 
led to such a formal technique, trying to identify the truly early sources of it. 
The historical course treated by us in pursuing this end, comprises two main 
steps, namely a first one having an introductory character which is devoted 
to a very brief outline of the main formal elements of scattering theory, and a 
second one just centered on such 1959 Tullio Regge paper. In achieving this, 

1 Often also said to be Watson-Sommerfeld-Regge transform. 
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a special attention has been paid to those many unavoidable moments con¬ 
cerning history of mathematics and mathematical-physics whose preeminent 
presence, along this historical route, cannot be evaded or neglected. But, let 
us say it immediately: we will not deal with the wide history of scattering 
theory, because our main aim just consists in directly arguing on the rising of 
1959 Regge work about the introduction of complex angular momenta and, 
above all, on the related mathematical methods there involved, pointing out 
just some historical aspects concerning these latter but that have nevertheless 
been quite neglected by the usual historical treatments of the subject. 

1. Scattering theory: a brief general overview 

Most of information concerning quantum systems comes from collision the¬ 
ory. Following^ (Rossetti 1985, Capitolo I), the early crucial moments of the 
history of quantum theory has seen basically involved celebrated collision ex¬ 
periments: to mention the main ones only, we recall the 1911 E. Rutherford 
experiments on elastic collision of a particles on atoms, the 1913-14 J. Franck 
and G. Hertz inelastic collisions of electrons on gas atoms and molecules, the 
photoelectric effect, the 1923 A.H. Compton scattering of photons on atomic 
electrons, the 1927 C. Davisson and L. Germer diffraction experiments and 
the 1934 E. Fermi experiments on neutron collisions; on the other hand, the 
modern elementary particle physics is also based on collision phenomena: for 
instance, the celebrated 1983 C. Rubbia experimental observation of interme¬ 
diate vector bosons essentially was a proton-antiproton annihilation process, 
which constituted one of the first experimental confirmations of the Stan¬ 
dard Model. In this paper, we are interested only in those historical aspects 
of collision theory which have led to the dawning of resonance dual models 
of string theory, with a particular attention to the many intersection points 
with history of mathematics. 

One of the chief notions of collision theory for putting together experi- Scattering 
mental data and theoretical issues, is that of scattering cross-section. The Tfuons 
differential scattering cross-section da = cr(fl)dfl is defined to be the number 
of scattered particles within the infinitesimal solid angle dQ (with f2 com¬ 
puted in polar coordinates 0 < 6 < tt and 0 < ip < 2 ti) per unit of time and 
per unit of incident beam of target particle. The total scattering cross-section 
is defined to be 



( 1 ) 


2 This textbook is one the main references that we have taken into account in drawing 
up this paper. Surely, (Rossetti, 1985) is the central reference considered in drawing up 
the Section 1. 
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which depends only upon the energy E of the incident beam. Following 
(Landau & Lifsits 1982, Chapter XVII) and (Rossetti 1985, Chapter II), in 
the center of mass reference frame, the motion of two particleqj with well- 
defined energy E and interacting by means of a central potential field V(r), 
may be reduced to the motion of a single particle with reduced mass m 
and relative energy E , which is ruled by the following well-known stationary 
Schrodinger equation 

h 2 

(2) - A -0(f) + V(r)ij}(r) = Eijjfr) 

2m 

which is the equation of motion of a particle with mass m moving in a central 
potential field V(r). If we set k 2 = (2 m/h 2 )E and U(r) = (2 m/h 2 )V(r), the 
equation (2) will assume the form 

(3) A'0(r) + [ k 2 — U(r)]^(f) = 0. 

We suppose U(r) —» 0 as f —> oo in such a rapid manner that we might 
speak of so afar regions where the interaction is negligible and where the 
above parameter k will be equal to the modulus of the wave vector k of the 
relative motion. 

We consider an incident beam directed along the positive z axis, as polar 
axis, oriented according to the related wave vector k. At infinite distances, 
the wave function 0 as (f) will satisfy the following asymptotic wave equation 

(4) A i/jas(r) + k 2 ip as (f) = 0 

whose solutions are either (progressive and regressive) plane waves of the form 
^asif) — Ae ±lk ' r and (divergent and convergent distorted) spherical waves of 
the form •0 as (f') = f( y 9,ip)e ±lk ' r /r, so that the wave function ip(r) describing 
a scattering process, at infinite distance, is the superposition of an incident 
plane wave (associated to the incident beam) and of a wave emerging from 
the scattering center. Therefore, we have the following asymptotic behavior 

(5) ip(r) ~ ip as (r) = A\e lkff +f(6,<p)- —) 

where A is a normalizing constant, f(6 , cp) is a function, called scattering 
amplitude and having a length dimension, which estimates the amplitude of 
the scattered distorted spherical wave with respect to the amplitude of the 
incident beam, 6 is the planar scattering angle of the emerging beam with re¬ 
spect to the z axis and k is the modulus of the wave vector k. The geometry 

3 With negligible spin-orbit interaction. 
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of the system implies that f(6 , <p) depends on the scattering angle 6 (colat¬ 
itude) but not on the anomaly (p . In the scattering amplitude are included 
the most important information on the scattering process, whose knowledge 
involves exact solutions to the original Schrodinger equation, hence the form 
of U(r ). Furthermore, as regard the differential scattering cross-section, it is 
possible to prove that 

(6) a(Q) = a{6,(p) = \f(0,<p)\ 2 
while, as regard the total scattering cross-section, we have 

( 7 ) a t (k)= [ \ne,<p)\ 2 dn, 

Jn 

in the geometry of central potentials having these formulas themselves but 
with f{6) for f(0,(p). Therefore, the scattering cross-section may be esti¬ 
mated by means of the scattering amplitude, and vice versa. 

In non-relativistic potential scattering theory, we have to consider a cen¬ 
tral short-range potential held U(r). Usually, the z axis will be the axis for 
the scattering center, coincident with the origin of the given polar coordinate 
system, parallel to the wave vector k — p/h of the incident beam. There¬ 
fore, we have a dynamical problem having cylindric symmetry with respect 
to the z axis, so that almost all the most important physical quantities will 
be independent from the anomaly tp (like f(6,ip) —> /(#)). In this case, un¬ 
der certain amplitude normalization conditions, equation (5) reduces to the 
following 

(8) #r) ~ + — 

r—>oo T 

which is a solution to the following Schrodinger equation 

(9) [A + (k 2 - U(r))]iP(r) = 0, 

and is the sum of two terms, an incident plane wave and an emerging spherical 
wave (scattered wave) that contains the scattering amplitude from which it 
is possible to deduce all the physical characteristics of the collision process. 
The equation (9), for central potentials, has a set of elementary solutions 
given by 

(10) il>im(r) = <p) 

which forms a complete set of solutions, so that every other solution to (9) 
will be a linear combination of them. In the geometry of the dynamical 
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problem that we have chosen above, only those spherical functions not de¬ 
pending on the anomaly ip may be considered, that is to say Y®(Q,<p) = 
\J(21 + \)/AtiPi (cos 0 ) , so that our wave function has a Legendre polynomial 
series expansion of the type 

OO 

(11) ip(f) = 0) = Cm(r)Pi (cos 9) 

1=0 

with pi(r) radial wave functions which have well-defined behaviors as r varies 
in [0, oo[. In particular, we are interested in the asymptotic behavior for 
sufficiently large values of r in which potential and centrifugal terms (due to 
finite potential barriers) are either negligible, given by 

(12) lMr) r^Tr An ( kr ~ l l + S ) 

where Ai is a suitable coefficient which goes to 1 as U{r) —> 0 and Si is an 
asymptotic phase displacement or shift of the given solution to the radial 
Schrodinger equation, related to the physical interaction and computed with 
respect to the free interaction solution. On the other hand, as U —> 0, the 
wave function defined by the asymptotic formula (9) should reduce to the 
free solution given by the incident plane wave, so that, taking into account 
(11), we will have 

OO 

(13) #r) —> V> ee (r) = e ikrcosd = £7(2/ + l)jf (fcr)^(cos 0) 

1=0 

where ji(kr ) are the Bessel spherical functions in the argument kr. In turn, 
(11) might be written as an expansion in partial waves each having a definite 
angular momentum /, of the typfEI 

OO 

(14) -0Cr) = + 1 )^( r )^( cos61 ) 

1=0 


where p>i(r) are solutions to the radial Schrodinger equation, regular at the 
origin and normalized in order to reduce them to Bessel functions ji(kr) as 
U 0. 

Inserting (12) into (14), we have 

OO A 

(15) ijj(r) ~ Y i\2l + l)-± sin (kr-£ + fcWcosfl) 

r->o ^' kr V 2 J 

1=0 


Phase 

shifts, 

partial 

ampli¬ 

tudes 


4 First expansions of this type were provided by W. Gordon in (Gordon 1928). 
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so obtaining an asymptotic partial wave development of the wave function 
which basically dependent on the phase displacements, or scattering phases, 
or phase shifts, Si, since it is possible to prove that also the constants Ai are 
function of the latter. In turn, Si are determined by the potential U (r) and 
vice versa: indeed, the knowledge of Si will allow the reconstruction of the 
interaction potential (scattering phase-function method). Therefore, these 
latter play a very fundamental role in potential scattering theor>0. They have 
real values for real values of the potential, whilst, for complex values of the 
potential, the phases have in general complex values as well. On the other 
hand, as pointed out above, the scattering cross-section is closely related 
with the scattering amplitude which, as now we shall see, is related with 
the asymptotic form of the wave function, hence, ultimately, with the phase 
displacements. Thus, measuring scattering cross-sections, we shall be able 
to get information about phase displacements, hence about the interaction 
potential. Indeed, following (Rossetti 1985, Capitolo II), it turns out to be 

-j OO 

m =5^ E< 2 '+bV" 1 - iwi ™ 5 ») = 

** 1=0 

1 

(16) =-V^(2Z + 0 e2lSl sin <5* • Pi(cos9) = 

rC 

1=0 

oo 

=E< 2 ' + l)aiPi(cos9) 

1=0 

having put, as a partial amplitude corresponding to the partial wave of an¬ 
gular momentum l 

(17) ai = ai(k) = ^(e 2l<5i - 1) = ^ e lSl sin <5,, 

5 Following (Landau & Lifsits 1982, Chapter XVII), the problem how to build up the 
form of scattering potential from the scattering phases supposed to be known, has a par¬ 
ticular relevance and has been solved by I.M. Gelfand, B.M. Levitan and V.A. Marcenko. 
In this regard, it turns out to be that, for the determination of U(r), it is enough to know 
5o(k) as function of the wave vector in the region comprised between k = 0 and oo as well 
as the coefficients a n in the following asymptotic expressions 

tpi ft aie^ 2m ^ ElUh)r /r 

of the radial wave functions <pi (of (11)) relative to those states corresponding to negative 
discrete energy levels, if these exist. The determination of U(r ) from these latter data 
involves the resolution of a certain linear integral equation (see (De Alfaro & Regge 1965)). 
Furthermore, another very interesting method to determine the scattering phase is the so- 
called variable phase method (see (Calogero 1967)). 
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so that, being a(6) = \f(9)\ 2 , we have the following expression for the differ¬ 
ential scattering cross-section 


( 18 ) 


a(6) 


1 | 


^(2 1 + l)e 2wSi sin <5; • Pi(cos6 ) 
1=0 


2 


whilst for the (elastic) total scattering cross-section we have 


/ I 4 oo oo 

a(6)d cos 6 = — ^(2/+l) sin 2 8 t = 47t ^(2/+l)|cp| 2 
1 1=0 1=0 

which, for 9 = 0, reduces to the optical theoren<Jtot(k) — (47r//c)9f/(0). 
From (17), it follows that, if the phase 5i = 5i(k ) assumes, for a certain energy 
valut0 of k, the value of 7r, then afk) = 0 and therefore the l- th partial wave 
does not concur to the scattering process (Rams auer- Towns end effect), whilst 
if the phase <5/ assumes, for a certain energy value of k, the value of tt/2, then 
the contribution of the /-th partial wave will be maximally exalted in the 
scattering process (resonance scattering). In the latter case, the scattering 
contribution by the l -th partial wave has the highest influence only when the 
phase 5i changes in a neighborhood of the energy E 0 corresponding to that 
value of k for which 8 [ assumes the value of 7t/2, in this case being possible 
to say that there exists an /-th partial wave resonance for E = Eq. In such a 
neighborhood, say X(Ef), an approximate expression for the /-th amplitude 
ai was provided by G. Breit and E.P. Wigner in 1946 (Breit- Wigner formula), 
and given by 


1 


where T = 2 /(d8i(E)/dE) e=e 0 - It provides a contribution to the total scat¬ 
tering cross-section given by 

(21) ai (E) = 4,(21 + 1)M 2 « ^(21 + 4 (g _ff +p/4 

whose graph in function of E shows a maximum value in the neighborhood 
X(Eq) having a width approximately equal to T. The appearance of an /-th 
wave function resonance peak physically corresponds to the production of a 
metastable (or resonant) state , localized in the interaction domain and with 
an equal angular momentum /, in which a particle may stay for a sufficiently 

6 Introduced by Melvin Lax in (Lax 1950a,b). 

7 Indeed E = h 2 k 2 /2m. 
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long time estimated by the mean lifetime given by r = h/T, in such a manner 
to exalt the corresponding scattering cross-section, this being the key essence 
of a scattering resonance phenomenon. The Born approximation for phases 
and scattering amplitudes respectively are 

poo poo 

(22) 5i = —k / U(r)[rji(kr)] 2 dr, f(0) — — / U(r) Sm ^ r r 2 dr. 

Jo Jo ( l r 

We are particularly interested in the complex case. When the poten¬ 
tial has real values, then the scattering phases are also real, but when the 
potential has complex values, then, in general, the scattering phases too are 
complex. In the case of central potentials, as regard the scattering amplitude 
f(9) given by (16), we consider the following formula 

OO 

— ]T(2i + l)(e 2 “'-l)P,(co S f>) = 

^ 1=0 
oo 

—5>( + i)(s,-i)a(co S f)), 

^ 1=0 

having put Si = e 2lSl . As regard the elastic differential cross-section, we 
have cr e i(9) = |/(0)| 2 , while, taking into account (19), for the elastic total 
cross-section, we have 



<7. 


_tot 
el 


(24) 


w=i3E( 2i + 1 )i 5 '- 1 i 2 = 

1=0 
oo 

^(2/ + 1)[1 - 2e~ 2m ‘ cos(eSM { ) + e" 4a,5i ]. 


7r 

Is 


1=0 


Instead, the inelastic differential scattering cross-section is given by 


(25) 


Oi 


_ OO oo 

W = p X> + !)[!- I^l 2 ] = p £( 2 < + !)[! - 


= -4a<5n 


1=0 


k 2 


1=0 


Therefore, in those scattering processes in which are present many alternative 
phenomena, like elastic and inelastic collisions, it is needed to consider the 
following expression for the total scattering cross-section 


27 T 


vtot(k) =<r*f (k) + (T in {k) = — Ypl + 1)[1 - 3ftSj] 


( 26 ) 


1=0 


2tt 


— 1 + 1)[1 - 2e~ 2QSl cos(29fW,)]. 


1=0 



The temporal evolution t 0 > t of a quantum system is defined through a 
unitary operator S(t,t 0 ) which acts on an initial state |s(t 0 )) producing the 
final state | s(t)) = S(t, to)\s(to))- In a scattering process, it is not necessary 
to know S(to,t) for arbitrary values of the time, but it is often enough to 
consider time intervals for to <C 0 and t > 0 when one supposes that the 
interaction is efficient only into a definite time interval centered around t = 0; 
in such a case, we may define a unitary operator S, called S matrix (or 
scattering matrix ) and introduced by J.A. WheeleiH in 1937, as a limit of 
unitary operators in the following fashion 

(27) S = lim S(t 0 ,t) 

£q —> — oo 
t—> oo 

where the asymptotic states to = —oo and t = oo refer to ideal moments in 
which the dynamical system is supposed to be into a non-interacting state. 
The S matrix (together the so-called T matrix, defined by S = I + iT ) is 
a formal tool which was mainly set up to study relativistic interactions of 
elementary subnuclear particles whose related forces were so little known in 
the form and nature to entail a quantum description through the Schrodinger 
and Dirac equations almost unusable. Instead, by means of such a tool, the 
scattering amplitudes could be deduced from general physical and mathemat¬ 
ical principles no matter by the knowledge of the involved forces. Further¬ 
more, the proliferation of the impetuous phenomenology of 1950s high-energy 
physics led to a so great amount of new and unexpected phenomena that the 
old quantum field theory turned out to be unable to provide a right theoret¬ 
ical basis to them, to make necessary appealing to general formal principles, 
above all unitarity and analyticity, to theoretically explain these new results. 
Every element of the S matrix, that is to say Sfi = (f\S\i), is the proba¬ 
bility amplitude to observe a particular final state |/) starting from a given 
initial state |i), that is to say, the probability amplitude associated to the 
dynamical process |i) —>■ |/). Other basic formal properties of the S matrix 
may be deduced from general physical principles, like relativistic invariants. 
For instance, as regard collision processes involving four scalar particles of 

8 Following (Mackey 1978, Section 21) and (Chew 1962, 1966, Chapter 1), the idea of 
using the S matrix (more properly, the S operator) to describe the ’’scattering” of particles 
by one another in quantum mechanics, was introduced by J.A. Wheeler in 1937 and again 
independently by W. Heisenberg in 1943 who however lost interest in such a tool because 
of formal problems related just to its analytic continuation that was required to give to 
the S matrix dynamical content. Later, C. Mpller published two seminal papers, in 1945 
and 1946, in which the mathematics underlying Heisenberg ideas were developed, so that 
sometimes S matrix is also called M0ller wave matrix. For a deeper historical analysis, see 
for instance (Chew 1962, 1966), (Cushing 1986) and (Mussardo 2007, Part IV, Chapter 
17, Appendix A). 
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the type 1 + 2 —y 3 + 4, the possible relativistic invariants are s = (p\ +p?) 2 , 
t = (pi — p 2 ) 2 and u = {p\ — P 2 ) 2 , which satisfy the following functional de¬ 
pendence relation s + t + u = Xu=i where m 2 = p 2 — — Pi ■ Moreover, 

another fundamental property of S matrix is the so-called crossing which re¬ 
lates the scattering amplitude of the initial process 1 + 2 —» 3 + 4 with those 
processes obtained by this replacing one incoming particle with momentum 
Pi of the former process with its antiparticle having four-momentum — pi, 
and vice versa. For example, in terms of T matrix, if T(s,t,u ) is the scat¬ 
tering amplitude of the process 1 + 2 —)■ 3 + 4, then it is also the scattering 
amplitude of the process 1 + 3—)-2 + 4orof the process 1 + 4 —>■ 2 + 3. 
Therefore, information about different collision processes may be suitably 
put into analytical relationship between them: for instance, one of these is 
given by that providing the presence of poles in correspondence to possible 
bound states or resonances. The scattering amplitude cannot have only poles 
because them generate only a real amplitude whilst the unitarity condition 
implies the existence of an imaginary part, say QT(s, Ujt), proportional to 
the total scattering cross-section, and that, in turn, contains the various 
cuts of the corresponding Riemann surface. Therefore, crossing and unitar¬ 
ity link together analytical properties of different collision processes, that 
is to say, different elements of the S (or T ) matrix, a notable fact, this, 
which was overstated in the 1960s till to think that such a formal schemqj 
could determine all the scattering amplitudes through well-determined spec¬ 
tral representations, amongst which certain dispersion relational, involving 
the various singularities of the analytic S (or T) matrix. 

Setting 


(28) 

we have 


|sj) = lim |s(t)), |s/) = lim |s(t)), 

£—>—00 t —^00 


(29) \s f ) = S\si), 

so that into the matrix S are included all the information on the scattering 
process by which the initial state |s,) temporally evolves to the final state 

9 Upon which relies the so-called bootstrap hypothesis according to which, roughly speak¬ 
ing, from a given set of bound states - named ’’particles” - describing physical states and 
forces through a well-defined dynamics formulated in such a manner that other states may 
be self-consistently generated without make reference to initial fundamental states. 

10 Following (Cini 1977), the dispersion relations had a preeminent role in high-energy 
physics, becoming the prevalent paradigm for strong interactions, within which frame¬ 
work a particular analytical property of the scattering amplitude related to non-relativistic 
Schrodinger equation was extended, by analogy, to the relativistic context just by Regge’s 
work. See (Cini 1977) for a deep historical account and basic philosophy of science con¬ 
siderations related to dispersion relations. 
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| Sf) due to the existence of an interaction between the particles constituting 
the given dynamical system involved into the scattering process. Now, if 
|Si) = |a), then the final state has the following expansion into a set of 
complete orthonormal states | b) 

(30) | S/ ) = f C at \b ) = [ <6| S/ >|6> = f <!>|S|a)|!>} = f SJ6), 

where the index a refers to the state evolution |a) —* |6), while S a b = (6|<S'|a), 
thanks to the unitarity of S, provides the probability amplitude, given by 
W ab = | S a b | 2 1 with which the dynamical system may run along scattering pro¬ 
cess from |a) to \b). In general, the final state | b) is, roughly speaking, char¬ 
acterized by a continuous variability in such a manner that it is more correct 
to formally consider a differential probability defined as dW ab = \(b\S\a)\ 2 db 
which, therefore, is the transition probability from the initial state |a) into 
one of the final states \b) having quantum numbers comprised between b and 
b + db. In the time-dependent interaction representation, the S matrix has 
the following representation 

j r°° 

-4 / m)dt 

(31) S = Te Tl d-oo 

said to be the Dyson’s formula, T being the time decreasing operator called 
chronological product and H' the perturbative part of the harniltonian op¬ 
erator. At the first order perturbative effects, this formula reduces to the 
following Born approximation for S 

j r°° 

(32) S^S B = I - - / H'(t)dt. 

nJ -oo 

In a scattering process, over each energy surface E , we have the following 
conservation energy relation 

(33) <&|S1a> = (b\S E \a)S(E b - E a ) = S E 5(E b - E a ) 

where S ba = (b\S E \a) is the so-called S matrix element over the energy surface 
E. Because of the continuous variability of the final states b, we have 

1 77 Born r\ 

(34) dW ba = —|<^|a>| 2 J(- -^\(b\H'\a)\ 2 Pf(b) 

which is said to be the first Fermi golden rule, where 

(35) Pf(b) = db/dEb 
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is the density of the final states. With greater precision, the two above 
mentioned quantum states, involved in the scattering process and considered 
satisfying the right normalization conditions, should be denoted as follows 
| a,E a , n a ) and \b,E b ,n b ), where a,b refer to the nature of the system and 
to the related characterizing quantum numbers, while n a = p a /p a — k a /k a 
and n b = Pb/Pb = kb I kb identify the related momentums whose directions are 
respectively identifiable by means of the related solid angles = (9 a ,ip a ) 
and Vlb = {9b, </?&)■ In particular, we have db = dE b dfl b , so that, in the 
hypothesis E = Eb = E a , the first right hand side of (34) becomes 

(36) dW ba = ^\{b,E,n b \S E \a,E,n a )\ 2 dn b 

which provides the transition’s probability per unit of time, from the initial 
state |a) to the final state \b) with wave vector having direction n b com¬ 
prised into the elementary solid angle dilb- Making use of the probability 
current density associated to a quantum particle, it is possible to prove that 
the differential scattering cross-section for the elastic and inelastic scattering 
processes a —> b is given by 

(37) a ba (il b )dil b = ^ = ^-\(b,E,n b \S E -I\a,E,n a )\ 2 dil b , 

\ja\ k a 

where j a is the probability current density associated with the particles of 
the incident beam. For elastic scattering processes, since a = b, we may write 

Ayr^ 

(38) a{ n) = o{0^) = —\( ( n'\S E -l\n)\ 2 

where = {6, tp) identifies the direction n' of the scattered wave vector with 
respect to the direction of the incident wave vector n. From (6) and (38), it 
follows that 

t 

(39) f{9, 9) = —e ia {n'\S E - I\n) 

K 

where e ta is a suitable phase factor. 

In a central potential problem, the angular momentum L is a constant of 
motion so that, said \lm) an orthonormal basis of eigenstates of the angular 
momentum operators L 2 and L z , the S matrix elements are diagonal and we 
may set 

(40) ( b , /W|S' E |a, Im ) = 5 w 5 mm 'S l ba 
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where the reduced S matrix elements S l ba (K) does not depend on m, seen the 
arbitrariness of the direction z with respect to which compute L z , but only 
on the two states a, b and on the energy E via the modulus k of the relative 
wave vector. The elements of the operator S E — I upon the energy surface 
E, given by 

(41) {b,n b \S E - I|a,ra 0 ), 
can be estimated as follows 

1 OO 

(42) {b,n b \S E - l\a,n a ) = — ^(21 + l)(S ba - S ba )Pi(cos9) 

71 1=0 

where 9 is the scattering angle between n a and n b , such a formula provid¬ 
ing the Legendre polynomial expansion of the matric elements (39), with 
coefficients which depend on the reduced S matrix elements S\ M . For elastic 
scattering processes, a = b so that, setting S l aa = Si, we have 

1 OO 

(43) (n'\S E - I|n) = — ^(2/ + 1)(5, - l)P,(cos0). 

n 1=0 

But, from (39) for central symmetry problems, we have 

2-7T 

(44) f(9) — —e ia {n , \S E — I|n) 
with aeRan undehned parameter, so that 

%ol 

(45) m= e —J2(2l + l)(S l -l)P l (co S 0) 

K 1=0 

which, compared with the following Legendre polynomial expansion 

1 OO 

(46) m = _^(2/+ l)(e 2 "' - l)P,(cos9), 

z=o 

implies e ia = — i and Si = e 2l&l , that is to say, the reduced S matrix elements 
are closely related to the corresponding phase displacements Si. Therefore, 
(39) reduces to 

Ojrn 

(47) f(9,tp) = - — {n'\S E -I\n) 

which holds too for all the elastic scattering processes, comprised the non¬ 
central ones, being 6 and ip the polar angles of n 1 with respect to the incident 
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direction n. In the case of central potential scattering, using (42) in (37), we 
obtain the following expression for the differential scattering cross-section 

1 OO 

(48) atafSlj) = jp| £(2i + !)(SL - Ufl(cosO)| 2 , 

1=0 

whilst for the total scattering cross-section, we have 


__ OO 

(49) 


both these latter formulas being valid for every scattering process a b. For 
elastic scattering processes (a = b), they respectively reduce to the following 


(50) 


and 



i p|£(2!+l)(S I -I)P,(co S <4)| 2 


OO 

(51) <rS*(fc; «->») = ££(2( + 1)|5, - I| 2 , 


whereas, for inelastic scattering processes (a ^ 6), we respectively have 


(52) 


and 


<%?(») = iF l£(2(+l)Sj a fi(co S f))| 

=0 


oo 

(53) <r“,(i; a -4 6) = p £(2i + 1)|SJJ 2 


Often, a sum over all the possible scattering channels a —> b is considered as 
follows 


(54) 


_ OO 

= £ <(k) = p £(2i + 1) £ |SjJ 2 

a^b 1=0 a^b 

__ oo 

^Da + Dd-W). 

1=0 
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Furthermore, in both elastic and inelastic scattering process case, we have 
Si = e 2lSl with \Si\ < 1, so that | S/\ = \e 2l5l \ = e~ 2 ^ Sl and A 5i > 0. As regard 
the total scattering cross-section, we have 

Q °° 

(55) o M (k) = a'f(k) + o‘Uk) = T7 £(2( + 1)(1 - «S,). 

1=0 

On the other hand, taking into account (46), we have 

OO 

(56) 3/(0) = —^(2( + l)(l-*S,), 

^ 1=0 

that, compared with (55), gives rise to the optical theorem, in its most general 
form 

(57) A/(0) = ^-G tot (k). 

In analogy with (47), which links together the elastic scattering amplitude 
with the elements of the S matrix, we may define a generalized scattering 
amplitude as follows 

(58) fba(n b , n a ) = f ba (^b) = —£-( bn b \S E - I| an a ) 

where fl b = O = ( 9 , (p) is the solid angle which identifies the direction of n b 
with respect to n a , this last chose parallel to the direction of the polar axis 
in such a manner that the differential scattering cross-section for the channel 
a —> b, is given by a ba (Q b ) = |/{, a (^b)| 2 - In the case of multichannel processes, 
for simplicity’s sake we consider a three channel process whose directions are 
n, n' and ft", by which we have 

(59) f{tf,n)-f*{n,f?) = ^ j f*{n"> n')f{n",n)dQ," 

that is the general condition to which undergoes the scattering amplitude for 
purely elastic processes, due to unitarity of S, and that, for central collisions, 
reduces to the optical theorem. Indeed, in such a case, f(n', n ) depends on 
cos# = n ■ nl only, and since we have too f(n,n') = it follows that 

(58) reduces to 

(60) $/(n,n') = i- J f*(n", n')f(n", n)dQ" 
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that is an extension of the optical theorem because, for n = n', we have 
9 = 0, so that (60) reduces further to 

(61) 3/(n,n) = J \f(n" ,n)\ 2 dVt" = 

that is to say, the optical theorem (57). It is also interesting to notice that 
(60) allows, at least in principle, to rebuild up the scattering amplitude when 
it is known its modulus, that is to say, it is known the differential cross- 
section, excepts the indeterminacy arising from the invariance of (60) with 
respect to the change f(9) —> —f*(9). 

Many fundamental results of scattering theory may be deduced without 
make reference to a specific form of the related interaction, only through the 
basic properties of the S matrix like, for instance, its unitarity from which 
we have information about the scattering amplitude. Further properties may 
be deduced from other analyticity properties of the elements of the S matrix, 
which are able to describe a collision process. In what follows, we shall refer 
to elastic collisions in a central held, where the S matrix elements are given 
by the functions Si = Sdk) - that is to say, S matrix elements computed 
over a surface of energy 11 ! k with values related to states of definite angular 
momentum l - depending on the scattering phases through 

(62) Si(k) = exp (2iSi (k)). 

For instance, for an elastic scattering process in a central held, if one considers 
(16) wrote in the form 

OO 

(63) f(k, 6) = YXU + l)ai{k)Pi(cos9), 

1=0 

where the partial amplitudes ai(k) are given by (17) and written as follows 

(64) a,(*) = d_(S,(A0 - l), 

then it follows that properties of f(k,9) are consequence of the properties 
either of Si(k) and of the Legendre polynomials P(cos9). The latter formulas 
(63) and (64) resolve, d’apres H. Faxen and J. Holtsmark works of 1927, the 
problem how to express the scattering amplitude by means of the scatter¬ 
ing phases. Moreover, from the orthogonality properties of these latter, an 
inversion of (63) gives rise to 

1 f 1 

(65) a>i(k) — — / f(k, 9) Pi (cos 9)cl cos 9 

2 J -1 

11 Recall that E = H 2 k 2 /2fi. 


Regge 

poles, 

Jost 

solutions 

and 

functions 
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hence, from (64), it also follows 


(66) Si(k) — 1 + ik J f(k,9)Pi(cos9)dcos9, 

so that there exists a certain equivalence between the informations given by 
the scattering amplitude f(k, 9) and those given by the set of functions Si(k). 
Therefore, it is of fundamental importance to study the analytical properties 
of the functions Si(k) considered as function either of l and k, though only 
for k G R + and l G No they have a direct physical meaning for scattering pro¬ 
cesses. Nevertheless, as we will see later, the consideration of these functions 
also for complex values of both these two variables, gives rise to interesting 
and unexpected results and interpretations. For instance, the knowledge of 
the properties of analyticity in the variable k of the functions Si(k), allows, 
inter alia, to show that the scattering amplitude verifies certain dispersion 
relations in the variable k, whereas the knowledge of the properties of analyt¬ 
icity in the variable l of the functions Si(k), provides a representation for the 
scattering amplitude which is useful to study the asymptotic properties for 
great transferred impulses, so giving rise to the so-called Regge poles. There¬ 
fore, for pursuing this, it is preliminarily need to define Si(k) for arbitrary 
values of / and k, and this may be done, in potential theory 12 l. through the 
following reduced radial Schrodinger equation 

(67) u '; + [e- l l±±l- U ( r )} Ul = o 

which, mathematically, has always a meaning for every complex value of k 
and l. The relation (64) has been obtained by taking into consideration a 
wave function with an asymptotic behavior given by (5) and developed in 
partial waves as follows 


( 68 ) 


VKO = Ih‘(a +1) 

1=0 


Xi(k,r) 

r 


Pi(cos9) 


where Xi(k,r) are the physical solutions to (67) which are zero in the origin 
of the given reference frame, and having the following asymptotic behavior 

(69) xi(k, r) ~ — — sin [At - Z J + Si(k)} 

r—too K Z 

12 Hence, within the non-relativistic quantum mechanics framework where an interac¬ 
tion’s potential is definable, so that the dynamical equation governing the collision process 
is the Schrodinger equation. 
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which, in absence of any interaction, reduces the following free solutions 


(70) 


X?(fc,r) = rj,(k,r) 


where ji(k,r ) are the Bessel spherical functions. The formula (69) implicitly 
defines the scattering phase Si(k), hence Si(k) via (62), through Xi(k,r) for 
each l,k. Thus, the analyticity properties of the functions Si(k) will be 
known once we know those of the functions Xi(k,r) as solutions to (67). For 
simplicity’s sake, we shall consider simple solutions to (67) as, for instance, 
those obtained under certain boundary (in the origin or at infinity) conditions 
which are very simple and which are dependent only on one of the two variable 
k and /, from which, therefore, to accordingly derive fundamental analyticity 
properties of Si(k) in the variable k (with l fixed at a given physical value) 
or in the variable l (with k fixed and belonging to M + ). If we set A = l +1/2, 
then (67) reduce to the following more symmetric form 

(71) u" + [k 2 - X ' ~ 1/4 - U(r)]u = 0 


whose solutions are of the form u = u( A, k,r ), or else u = ui(k,r ) when we 
wish to refer to physical values of l. We search for solutions to (71) through 
suitable hypotheses on the interaction potential U(r) which, nevertheless, 
allow us to have an enough degree of generality to leave aside from the 
particular type of potential so chosen. Such hypotheses regard the behavior 
of the potential either in the origin of the reference frame and at infinity. To 
be precise, we consider short-range potentials U (r) such that 


(72) 


r\U(r)\dr —M(a) < oo Va G 


\U(r)\dr —N(b) < oo Wb e 


The hypothesis (72) i implies that the point r — 0 is a singular point of 
Fuchsian type for the equation (71), with characteristic exponents A + 1/2 = 
l + 1 and —A + 1/2 = —l, so that surely there exist two solutions to (71) 
which behave like to r ±A+1 / 2 as r —> 0. Instead, the hypothesis (72)2 implies 
that (71) has solutions which have an asymptotic behavior similar to the 
one of the solutions to the equation u" + k 2 u = 0, so that surely there exist 
solutions to (71) which behave like to exp (±ik ■ r) as r —> oo. Therefore, 
through such asymptotic behaviors, namely r ±A+1 / 2 as r —> 0 and e ±lk ' r as 
r —> oo, we are able to find solutions to (71). 
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The solution u = 0(A, k,r) ~ r A+1 / 2 , with the boundary condition 

r —>0 

(73) limr -A-1/,2 0(A, k, r) — 1, 

i— 

is said to be the regular solution to (71), and, as such, it turns out to be given 
by <j)(\,k,r) = (pi(k,r) = Ci(k)xi(k, r). Likewise, the function 0(A ,k,r), in 
absence of interaction, must reduce to the free solution 0°(£;,r), that is to 
say, to the solution to the free equation 

(74) u » + (e - OLzllly = o 

which has the same asymptotic behavior of 0(A, &,r) in r = 0. For l e No, 
we have 

(21 4- IV 

(75) 0°(M = ^ " rMk,r) 

where ji(k, r) are the spherical Bessel functions. Therefore, to may carry out 
an extension to every value of /, we should consider the functional gener¬ 
alization of the semi-factorial, hence the Euler gamma function, so that we 
have 

(76) 0°(A, k, r) = 2 a £T a T(1 + A )r 1/2 J x (kr) 

where J\{kr ) = Ji+iMkr) are the spherical Bessel functions of the first kind, 
defined for arbitrary values of l. To find analytical properties of </>(A, k, r) as a 
solution to (71), we may construct a suitable integral equation of the Volterra 
type to which such a (j) must satisfy under the boundary conditions (72), and 
this, in turn, may be accomplished by means of an auxiliary equation that 
takes into account the same boundary conditions (72), equivalent to 

(77) 0 ~ r A+1/2 , 0' ~ (A + l/2)r A-1/2 , 

i—>o i—>o 

and that has known solutions having the same behavior, in the origin of the 
reference frame, of the searched solutions to (71), that is to say, having a 
behavior as r ±A+1//2 . As a possible auxiliary equation to which 0 ought to be 
satisfied as r —» 0, we choose the following one 

(7S) u" - A2 ~ 1/4 u = 0, 

which has the following pair of fundamental solutions 
(79) u+ = r A+1/2 , w_ = r" A+1/2 . 
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It is possible to prove that such an integral equation is as follows 


0(A, k,r ) = u + (r)+ 


(80) 


1 

2A 


[. k 2 — U(x)][u + (x)u-(r) — u-(x)u + (r)](j)( A, k,x)dx, 


so that, the integral equation of the Volterra type to which the regular solu¬ 
tion 0(A, k,r) must satisfy is 


<f>( A, k,r) = r A+1 / 2 + 

^ ^ + ^ ~~ ^ (“) }^(A,A;,a?)da? 

whose formal solution is 

OO 

(82) 0(A ,k,r) = ^T<l> n (r) 

n= 0 

with 0o ( r ) = r A+1 / 2 and 

(83) (j> n (r) = ^ I [k 2 - C/(x)]x/rx| ^ Q }0 n _i(x)ch, Vn e N, 

which provides too an effective solution to (80) when it is uniformly con¬ 
vergent. To study the related convergence, we distinguish between the two 
cases 9ftA > 0 and 9ftA < 0. Putting A — n + ia, in the case 9ftA — n > 0, 
which comprises the physical case of angular momentum values (which are 
semi-integers), it is possible to prove that 


(84) 


with P{r) 


\M r )\ < 


r (i+i/2 ^p^ r y n 

|A| n n! 


Vn e N, 



U(x)\dx < (1 /2)|A: 2 |r 2 + M(r), so that we have 


(85) |0(A, k,r)\ < £> n (r)| < r M+1/2 exp(P(r)/|A|), 

71=0 


a condition which guarantees convergence. Now, a celebrated theorem, due to 
H.J. Poincare, which roughly states that a solution to a differential equation, 
like (71), whose coefficients are entire functions of a certain parameter, is also 
an entire function of the same parameter when defined through boundary 
conditions which are independent from this parameter, may be applied to 
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(81). Indeed, (81) is, in every finite domain, a uniformly convergent series of 
polynomial in the variable k 2 , so that it is an entire function of k 2 as well, in 
agreement with the just above mentioned Poincare theorem. Moreover, it is 
possible to prove that, at least in the region 9ftA > 0, each term of (81) is an 
analytic function in the variable A, so that also 0(A, k, r ) is analytic there. 
The analytic continuation of 0(A, k, r) for negative values of 9ftA entails more 
detailed information about the potential U(r); in particular, if one supposes 
that 

(86) U(r) ~ cr~ 2+e , e > 0, 

r—>0 

then it is possible to prove that 0(A, k, r ) is analytic in A in the region 
9ftA > —e/2, so that there exists, at least, a region of the A plane, i.e. the 
strip |!RA| < e/2, in which the functions 0(A ,k,r) and <f>(—\,k,r) are con¬ 
temporaneously defined and constitute a fundamental pair of solutions to 
(71) for A / 0. Therefore, for each A arbitrarily fixed, </>(A ,k,r) is an entire 
function of k 2 , at least for 9ftA > 0, or else for 9ftA > 0 if (86) holds, as well 
as it is an even function of k, i.e. <f>(\,—k,r) = cj)(\,k,r)‘ furthermore, for 
each k, </>(A, k,r) is an analytic function of A at least for 9ftA > 0 or else for 
9ftA > —e/2 if (86) holds, for each finite value of k arbitrarily fixed. Instead, 
as concerns the behavior of </>(A, k,r) as k —» oo, it is possible to prove that, 
for any direction by which k —> oo, 0(A ,k,r) reduces to the following free 
solution 

(87) k, r ) = + c 2 |^) (1 + 0(l/k)) 

where Ci,C 2 are independent from k. Finally, we notice that 0(A ,k,r) is real 
when k and r are real (under the obvious hypothesis that U (r) is also real) 
because of the boundary conditions (73), so that, due to the H.A. Schwarz 
reflection principle, the following Hermiticity condition holds cp*(X,k,r ) = 
</)(A *,k*, r) which, for physical values of l = A — 1/2, implies </>*( k,r ) = 
4>i{k*, r ). 

We have seen that the S matrix elements, as well as the scattering phases, 
are determined by the behavior of the wave function as r oo, so that 
a primary role is played by boundary conditions at infinity, like (72); in 
particular, the condition (72) 2 says that U(r) goes to zero faster than 1 jr as 
r —> oo, and guarantees the existence of two independent solutions to (71), 
which, as r —y oo, behave as exp (±ik- r). Let /(A, k, r ) be the exact solution 
to (71) with the following boundary condition 

(88) lim e*^' r /(A, k, r) — 1, 

r—>oo 
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so that 


(89) f(\k,r) ~ <r*\ 

r—> oo 

which is said to be the Jost solution. Moreover, due to the inversion symme¬ 
try invariance of (71) with respect to the variable k, it follows that /(A, —k, r) 
too is a solution to (71), with 

(90) /(A, —k,r) ~ e* fcV 

r—> oo 

Therefore, for every k 0, we shall refer to the pair of Jost solutions 
/(A, ±k, r) as two independent fundamental solutions to (71). The free Jost 
solution is given by 


(91) 


/o(A, fc, r) 



/Q\ 

where 77/ ; is the second-kind Hankel function. Since (88) does not depend 
on A, the above mentioned Poincare theorem states that, for each k fixed, 
f(X,k,r), where is defined, is an entire function of A 2 , hence everywhere 
analytic over the finite plane of A. As an entire function of A 2 , f(X,k,r) 
is even in the variable A, that is to say, f(—X,k,r) = f(X,k,r). Indeed, if 
one proceeds as made above, the analytical properties of f(X,k,r) may be 
deduced transforming the differential equation to which it satisfies, i.e. 


(92) 


/" + 



A 2 - l/4-i 
72 


0 


with the boundary condition (88), into an appropriate integral equation of 
the Volterra type, through the following auxiliary equation 


(93) 


u" + k 2 u = 0 


and its solutions 
(94) u + = e 

Such an integral equation is 


ik-f 


u_ = e 


—ik-f 


/(A, k,r) =e~**+ 
(95) 1 k-oo 

+ k.l 


r A 2 — 1/4 

sin[fc(x — r)] | U{x) 4--- 


x* 


/(A, k, x)dx 
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from which it follows that /(A, k,r ), for each k ^ 0, is an entire function of 
A 2 reasoning as made above by means of a series expansion whose general 
n-th term is a polynomial of degree n in the variable A 2 , at least in the region 
fsk < 0. Likewise, it is possible to prove that, for any direction by which 
\k\ —* oo, we have 


(96) f(X,k,r) ~ /o(A, k,r). 

|fc|—>-oo 


The analyticity region in k may be extended if, for instance, one considers 
potentials U(r) such that 


/•OO 

(97) / rU(r)e ur dr < oo for v < m, 

Jo 

where m is an arbitrary positive number, then it is possible to prove that 
f(X,k,r ) is analytic in k in the region 3 k < m/2, except k — 0, which is 
a polydromy point having kinematical nature because it is also shared by 
/o(A, k, r), besides to become a pole of order l for physical angular momenta 
(l E No). Taking into account this, it is usual to cut the complex pland^l 
of the variable k along the imaginary axis between the origin and the point 
im/ 2, so giving rise to the so-called kinematical cut (because independent 
from the potential). Therefore, under condition (97), the function f(X,k,r) 
is analytically defined in the region Qk < m/2, hence /(A, —k, r) in the region 
fsk > —m/2. Furthermore, for Yukawian potentials of the type 


(98) 


U(r) 



e Air c(/r)<ip, 


it is possible to prove that Jost solution f(X,k,r) is, for every A arbitrarily 
fixed, everywhere analytic in the variable k, except in the branch point in 
k — 0 and in a cut in the positive imaginary semi-axis made between im /2 
and ioo, and said to be the dynamical ctiO Therefore, to sum up, the 

13 Following (Eden et al., 1966, Preface), one of the most remarkable discoveries in 
elementary particle physics has been that of the existence of the complex plane. From the 
early days of dispersion relations, the analytic approach to the subject has proved to be 
one of the most useful tools. 

14 For certain types of potential, like the Yukawian one (see later), such a dynamical 
cut may degenerate into a series of poles which are usually said to be false poles of the 
S matrix, and that, in general, are poles of fi(k) in the region > 0 where, usually, 
fi(k) is not analytic. Nevertheless, we are interested in the analytical properties of the 
scattering amplitude /(fc, d) that, in potential theory, does not have such a dynamical cut 
in the region > 0, so that every pole of f(k, 6) is there always associated with a bound 
state having a definite value of l (see later). 
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Jost solution is defined by the boundary condition (88), it is a holomorphic 
function in the finite plane of A (at least, in the strip |J?A| < e/2), is an even 
function in A and an entire function in A 2 . As regard its dependence on k, 
such a function is analytic in the region Ah < 0 and continuous in the real 
axis except a polydromy point in the origin which becomes a pole of order l 
for physical values of A = l +1/2. It is analytic in the region Ah < m/2 if the 
potential verifies (98); in particular, if m = 0, it turns out to be continuous 
in the real axis (except h = 0), while if m — oo, the unique finite singularity 
is the point h = 0. For Yukawian potentials, it is analytic in the whole of 
the plane h cut between the origin and ioo; for physical angular momentum 
values, such a cut reduces to the dynamical one between im/2 and ioo, plus 
a pole in h = 0. Moreover, thanks to (89), /(A, h, r) is, for real A and pure 
imaginary h = ir, to a real function of A and r, so that, due to the Schwarz 
reflection principle, it has the following Hermiticity property 

(99) r(A,h,r) = /(A*,-h*,r) 

that, for physical values of the angular momentum, reads as 

(100) ff(k,r) — /)(—h*,r). 

Therefore, in the above mentioned regions, the four fundamental solutions 
to (71), that is to say, the two regular solutions 0(±A ,k,r) (which are also 
functionally independent for A / 0) and the two Jost solutions /(A, ±k,r) 
(which are also functionally independent for fc / 0), are all well-defined, so 
that the ones are expressible as linear combination of the remaining others. 
We shall consider the regular </>(A, k, r ) as a linear combination of the Jost 
solutions /(A, k,r ) because, through it, we shall be able to have a link with 
the S matrix elements. To be precise, we consider the following combination 

(101) (j)(X,k,r) = 2^[/j(A, k)f(X, —k, r) - fj(X,-k)f(X,k,r)] 


where 

( 102 ) 


fi(k ) = fj(X,k) =w[f(X,k,r),<f>(X,k,r)] = 

—/(A, k, r)4>'(X, k, r) - /'(A, k, r)<f>( A, k, r ) 


is the Wronskian between /(A, k,r) and 0(A, k,r), which define a new func¬ 
tion said to be the Jost function. From (101), it is possible to prove that the 
asymptotic behavior of the regular solution is given by 


(103) 


0(A, k, r 


~ tyh\f J ^X) k)e lk ' r — fj(X, —k)e~ lk ' r ]. 

r —>oo Ilk 
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Moreover, by means of the boundary behaviors of the regular solutions before 
considered, it will be possible to attain the following expressions for such a 
Jost functional, namely 

/(A, k) =2A lim r x ~ 1/2 f( A, k, r) (3JA > 0), 

r—>0 

(104) 

JiXJ.-j =2/7- lim Cil; < 0), 

r—> oo 


thanks to which it will be possible to infer the properties of /(A, k ) from the 
boundary behaviors and the properties of the regular solutions. For instance, 
the analyticity properties of the Jost function, in the variable A, are the same 
of those of the regular solution </>(A, k,r ), whereas the analyticity properties 
of the Jost function, in the variable k , are the same of those of the Jost so¬ 
lutions f(X,k,r). In particular, we have the following Hermiticity condition 
/*(A, k) = f(X*,—k*) (which holds for real values of U(r) and r) that, for 
physical values of the angular momentum, implies f*(k) = fi(—k*). In ab¬ 
sence of interaction, from the free expressions for regular and Jost solutions, 
we deduce the following free Jost function 


(105) 


/o(A, k) =v/27vr2 A r(l + X)e~ i{x - 1/2) % k~ x+1/2 = 
=2A lim r A-1 / 2 / 0 (A, k, r ). 

r—»0 


We have above considered the regular solution 0(A, /, k) that, for physical 
values of A — l + 1/2, that is to say, for semi-integer values of A, it vanishes 
in the origin of the reference frame as r l+1 , so that it turns out to be propor¬ 
tional to the physical solution xi(k,r), that is to say cp(X,k,r) = (pi(k,r ) = 
Ci{k)xi{k,r), where Xi is the so-called physical solution, in the origin, to 
(67), defined by (69) that here we re-write in the following form 


(106) 


Xi(k,r) ~ 

r—>■ oo 


e lSl 
2 ik 




Ak-r 


_ e -i(Si-ln/2) 



so that we have the following asymptotic behavior 


(107) 


4>i{k,r) ~ 

r^-oo 


Ti{k) 

2 ik 


■ e i(5i-ln/2) e ik-r _ e ~i(5i-lir/2) e ~ik-r-\ 


15 Often, the expression /(0,fc) is also considered as a Jost function (see (Umezawa & 
Vitiello 1985, Chapter 7, Section 7.8)). 
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having put Ti(k) = Ci(k)e lSl . Considering complex values of / = A — 1/2 and 
of k, we re-write (107) as follows 


(108) 


M k ,r) 


r(\,k) 


rsj 

r—Hx 2 ik 


gi(i5(A,fc)—(A—l/2)7r/2) _ 


_|_ e -i(<5(A,fc)-(A-l/2)ir/2) g-ifc-rl 


and comparing this asymptotic behavior with the one given by (103), we 
have the following identities 

f(X,k) = r(A,A;)e i(5(A ’ fc) - (A - 1/2)7r/2)) , 

(109) 

/(A, -k) = r(A, /c)e _i ( 5 ( A ’ fe ) _ ( A_1 / 2 ) 7r / 2 )), 

which can be considered as a parametrization of the Jost functions /(A, ±/c). 
Making the ratio between the two formulas of (109), we see that the function 
S(X,k) = e 2i(S ( A ’ fc ) is linked to the Jost functions /(A, ±k) by the simple 
following relation 

(110) S( A, k) = 

J (A, —k) 

that, for physical values of /, reduces to 

(111) s,(k) = (-1 

M-k) 

which represents the analytical continuation, to complex values of k and /, 
of the elements Si(k) of the S matrix. Therefore, the S matrix elements 
S( A, k), as functions of the complex variables A and k, are determined by the 
Jost functions /(A, ±k). In particular, we are concerned with the properties 
of S( A, k) — S(l + 1/2, k) = Si(k ) for fixed physical values of / (and arbitrary 
k), and for fixed arbitrary values of k (and arbitrary A), as well as in their 
implications for scattering amplitudes. As concern some first analyticity 
properties of S( A, k), we may immediately say that S( A, k) is, for each hxed 
A with JJA > 0, a meromorphic function of k at least in the region \5sk\ < mj 2 
(in the hypothesis that U(r) goes to zero as r —>■ oo at least like e~ mr ) except 
the point k — 0, while S( A, k) is, for each fixed k (in particular, for physical 
values of k), a meromorphic function of A at least in the region JiA > 0. 
Considering Si(k) as a function of the complex variable k with a fixed / equal 
to one of its possible physical value, we have that (111) holds, from which it 
follows that, in the k plane, there may be two types of singularities: the first 
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type corresponding to zeros of fi(—k), the second type corresponding to the 
poles of fi(k) in the region 5sk > 0 (or 5sk > m/2). But what we want to 
highlight in this regard is the main fact that the zeros of the Jost function, 
hence some poles of Si(k), are determined by the nature of the dynamical 
system involved in the collision process, hence by its potential. For bound 
states, having a certain angular momentum l and defined by a radial wave 
function regular in the origin of the reference frame, say (f>i(k, r), we have that 
the latter is a square-integrable function (condition characterizing a bound 
state that, physically, corresponds to a particle which is relegated near the 
scattering center) only for k not real but purely imaginary, that is to say, for 
k = in with k e M + (take into account the symmetry of (f>i(k, r ) with respect 
to k —y —k). From (103), it follows that 

( 112 ) <M*«,r) ~ -^7 [fi (*«)<e“ Kr -//(“*«)' e * r ] 

r-j-oo ZK 

so in order to have a bound state for k = in with energy E = —K 2 k 2 /2/j,, it 
is necessary that the coefficient of e Kr is zero, that is 

(113) = 0, 
so that 

(114) <&(*«, r) ~ 

which is normalizable, whence it describes a bound state. Thus, (113) pro¬ 
vides a necessary and sufficient condition for the existence of a bound state 
having angular momentum l and energy E = —h 2 K 2 /2/i 1 that is, a bound 
state corresponds to a zero of the Jost function in the half-plane 5sk < 0 on 
the negative imaginary semi-axis Rk = 0. In such a case, the wave func¬ 
tion <pi(iK,r ) describing such a state is a real function because, taking into 
account the parity symmetry k — y —k, we have 

(115) (f>*(iK,r) = r) = <f>i(iK,r). 

Therefore, as regard the S matrix elements, the bound states manifest them¬ 
selves both as poles of S)(/e) in 5sk > 0, due to the zeros of fi(—k), and as 
zeros (placed symmetrically with respect to the poles) of fi/k) in 5sk < 0 
and placed into the semi-axis Rk = 0. It is also possible to prove that, for 
bound states, the latter are simple, so that the corresponding poles of Si{k) 
are simple too. If one denotes with fi(k) the first derivative of fi(k) with 
respect to k, then it is possible to prove that, for bound states, the following 
relation holds 

r°° i 

(116) J [</>i(ko,r)] 2 dr = -^fi(k 0 )fi(-k 0 ), 
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while, an expression for the residues of Si(k), with respect to the variable 
k = in (k > 0), into the poles corresponding to bound states, is as follows 


(117) 


Ri{ik) = -i(- 1 )' 


Nf{n) 


(jyf(in, r)dr 


where Ni(k ) = 2k is the normalization constant of the asymptotic 

wave function given by (114). It follows that such a residue is purely imag¬ 
inary, while, if one uses the variable E = h 2 k 2 /2/j, instead of k, in regard 
to a bound state corresponding to k = in, the residue itself computed with 
respect to E = E K = — h 2 K 2 /2fi is given by 

j L- Pj 2 

(118) r(E K ) = lim (E - E K )Si(k) = - Ri[ik) 

E^rE n /J, 

from which it turns out to be real. 

Besides the above mentioned poles in the half-plane c 2sk > 0, with respect 
to the parity symmetry k —> —k, Si(k) has further poles also in the half-plane 
5s k < 0 , corresponding to the possible zeros of fi(—k) in such a half-plane, 
where, in general, it is not analytic. In such poles, the wave functions <pi(k, r ) 
have an explosive behavior and, therefore, do not represent bound states. In 
the region £ sk < 0, we may have poles either with 97 fc ^ 0 and 37/c = 0. In 
the first case, if k = fco is a pole of Si(k) in the region 5s k < 0 , the case 
37fc 0 7 ^ 0 (which cannot to take place for bound states due to the related 
square-summable finiteness condition) implies the existence of pairs of poles 
symmetrically placed with respect to the axis 37fc = 0. The poles of Si(k) 
as zeros of fi(—k) = 0 in the region 5sk < 0 with 3 R,k ^ 0, are usually 
called resonance, while the poles of Si(k) as zeros of fi(—k) = 0 in the same 
region but with 3 R.k = 0 , are said to be virtual states (or antibound states 
or redundant zero J^ll . For a resonance, that is to say, for a pole of Si(k) of 
the type k = h + ib with b < 0 , we never have b = = 0 since, if it were 

fi(—k ) = 0 for k G M, due to the Hermiticity condition f*(k) = fi(—k *), we 
would have too f*(—k ) = fi(k) = 0 , so that the wave function <pi(k,r ) would 
turn out to be zero as well. We give the name resonance to such a type of 
pole because, when b is small, its presence is manifested in the scattering 
amplitude through the appearance, in the Z-th wave function component, of 
a resonance, centered around an energy value E = E 0 = h 2 h 2 /2[i and having 
a width T proportional to | 6 |. This last reasoning holds too for all the poles 
of Si(k) in 7y k < 0 with 3 R.k ^ 0, for every physical value of l. Instead, an 

16 The redundant zeros were first discussed by S.T. Ma in (Ma 1946; 1947a,b). 
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antibound state meant as a zero of fi(—k) in the region 5sk < 0 with Rk = 0, 
say k = ib with b < 0, may manifest with an exaltation of the scattering 
amplitude at the threshold, if |6| is very small, even if this last phenomenon 
is also due to the presence of a bound state. For what follows, in studying 
analytical properties of Si(k) for physical values of /, it is better to choose 
E = h 2 k 2 /2/i as a fundamental variable rather than 

(119) k = VE^2fi/h 2 . 

Therefore, to get full variability of k, E must variate over all the two sheeted 
Riemann surface of \/~E , made by two sheets with branch point at the origin 
except for a cut on both sheets extending along Rk < 0. In such a case, it 
is a common practice to cut the plane of the complex variable E along the 
positive real semi-axis, calling physical sheet (or first sheet ) of such a Riemann 
surface the one corresponding to the positive determination of \/~E (hence 
that corresponding to 5sk > 0), and non-physical sheet (or second sheet ) the 
one corresponding to the negative determination of \[E (hence corresponding 
to 5sk < 0). Thus, for Si = Si(E), we have that bound states are poles of 
Si(E) placed on the negative real semi-axis of the physical sheet (being, for 
bound states, 5sk > 0 and Rk < 0); the resonances are poles of Si(E) which 
lie on the non-physical sheet (being 5sk < 0) and the nearer they are to the 
cut, the stronger their effects will be on the scattering cross-section; and the 
antibound states are poles of Si(E ) which lie on the negative real semi-axis 
of the non-physical sheet (being 5sk < 0 and Rk < 0) and the nearer they are 
to the origin of the reference frame, the stronger their effect will be on the 
scattering cross-section in threshold; in this latter case, moreover, we have 
trajectories which behave for a while as bound states but finally break-up 
into two particles moving apart (see (Mackey 1978, Section 21)). Therefore, 
it is meaningful to speak of either close polar singularities with respect to 
the physical region E > 0, represented by the cut, and far polar singularities, 
the first ones being those that mainly contribute to determine the form of 
the scattering cross-section. Furthermore, if the pole of S/fk) is at the point 
k = a+ib, then it represents our unstable bound state of energy a and angular 
momentum (or spin) /; the larger b is, the shorter the particle lifetime will 
be, and when the lifetime is so short that the bound state decomposes before 
it can travel a measurable distance, then we have a resonance. Of course, 
there is no sharp distinction between resonances and unstable particles. It is 
interesting that resonances and unstable particles not only have energies and 
angular momenta but behave like stable particles in other respects. One can 
even discuss their interactions with one another (see (Mackey 1978, Section 
21 )), 
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Now, we consider the scattering amplitude f(k, 9) (as, for instance, given 
by (63)) as a function of the energy E = h 2 k 2 /2/i and of the scattering angle 
9 , hence denoted by f(E, 9). We may consider E as a complex variable. For 
simplicity’s sake, we consider natural units, so that h — 2p = 1 , whence 
k = \/~E. Therefore, considering f(E,9 ) as a function of E at each fixed 
value of 9 , from (63) we have the following partial wave development 

OO 

(120) f(E,6) = ^(2Z + l)a,(fc)P,(cos0) 

o 

where 

(121) a,(k) = T( Sl (k) - 1) 

and k = \J~E. Therefore, the singularities of the partial amplitudes a; are 
the same of Si(k), and since it is possible to prove that the dynamical cut of 
Si(k) does not appear in the scattering amplitude, we have that the unique 
singularities of f(E,9), in the variable E, are poles corresponding to the 
possible bound states, resonances and virtual (or antibound) states, besides 
the possible polydromy points of \/~E, i.e., origin and point at infinity. We 
restrict to consider f(E,9) over the physical sheet of \/~E , where the unique 
possible singularities, for what has been said above, are the points of the cut 
along the positive real semi-axis as well as possible bound states having a 
definite angular momentum l and placed into the non-positive real semi-axis. 
Moreover, the singularities of the scattering amplitude f(E,9) on the phys¬ 
ical sheet coincide with the spectrum of the Hamiltonian of the dynamical 
system constituted by the interacting particles, including the allowed bound 
states as poles of f(E, 6) for §IE < 0, which are in a finite number for every 
potential decreasing faster than 1/r 2 (hence belonging to the discrete spec¬ 
trum), and all the states of the cut of f(E, 6) which correspond to all values 
of ?R,E > 0 (hence belonging to the continuous spectrum). Then, as E —> oo 
by means of any direction along the physical sheet, it is possible to prove 
that f(E,9) f B {9) in a uniform manner, being fs{9) the Born scat- 

E—> oo 

tering amplitude obtained considering the Born approximation (22) 2 which 
holds when high energy values are called into question. To this point, it is 
important to show that the scattering amplitude f(E , 9) undergoes a certain 
energy dispersion relation computed with respect to a fixed value of the angle 
9. To get it, we consider the Cauchy integral representation of f(E , 9) — fs(9) 
extended along a closed path of the E' plane this last being cut along the 
dtE' > 0, having poles along $IE' < 0, surrounding the point E but no one 
of the singularities of the integrating function. Thanks to Cauchy theorem, 
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this integration path, say T, may be homotopically deformed into a bigger 
circle, with radius R —> oo, closely embracing the cut dtE' > 0 plus small 
neighborhoods each of which centered around every pole of dtE' < 0. As 
we have said before, f(E,9) — f B (9) uniformly tends to zero, so that the 
contribution of the integral 

( m) 

along this circle, goes to zero as well, so that we may write 


(123) 


f(E,9) =f B {0)+ 

1 f°° f(E' + ie,e)-f(E' 

2ni J 0 E' — E 


ie, 9) 


dE’ + Y .; 


Vn(0) 

E- En 


the hrst contribution being extended to that part of the path T embrac¬ 
ing dtE' > 0 along a narrow strip of width e —» 0, whilst the second one 
comprehends each contribution given by the neighborhood of every pole of 
$tE' < 0, where y n {9) = Res {f(E,9)} E=En , for each n, is the residue of the 
scattering amplitude in the pole E = E n (which represents a bound state 
having a definite angular momentum /) and the summation is over all the 
allowed bound states. For 9 e R, we have f*(E,9) = f(E*,9 ), whence 
f(E — ie, 9) = f*(E + ie, 9), so that (123) reduces to 


(124) 


f(E,e) = f B (e) + - 

7r 



mE’,8) 

E 1 -E 


dE' + Y 


n 


yn{9) 
E — E n 


where A f(E',9 ) is evaluated upon the upper boundary of the cut. The 
dispersion relation (124) is due the classical form of the Khuri’s dispersion 
relation (see (Khuri 1957)). This last expression is nothing but a dispersion 
relation for the scattering amplitude which provides all the possible singu¬ 
larities and characteristics: namely, the cut along the positive real semi-axis 
and the discontinuities of f(E, 9) in it (which are 2 i$sf(E, 9)), as well as the 
the poles corresponding to bound states for E = E n {'RE n < 0) and the 
related residues y n (9)', furthermore, it provides the behavior of f(E,9) as 
E oo, being f(E,9) —>■ f B (9). As regard, then, the residues y n {9 ), these 
can be evaluated through the corresponding functions Si(k) as follows 

(125) y n (9) = (2 1 + D ^ftU cos*) e R 

2 y Fj n 

where ri(E n ) is the residue of Si(k) in E = E n , basically given by (117) and 
(118). Nevertheless, the (124) not always turns out to be useful because it 
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might be quite difficult to have information on $ sf(E , 9), so that, through the 
optical theorem (57) in the energy E with respect to 6 = 0, another formula 
is often used instead of (124), namely the following one 


(126) 


s/(£,0) = / b (0) + 21 



E' — E 


dE' + Y. 


n 


Vn( 0 ) 

E - E n 


which is said to be the forward dispersion relation of the scattering ampli¬ 
tude, connecting the real part of f(E, 0) with an integral over its imaginary 
part as well as with the position and residues of the poles. It was provided 
by D. Wong and N.N. Khuri in 1957. 

We reconsider S(X,k), as given by (110), as a function of the complex Re gge 
variable A, for any k G M + arbitrarily fixed. For k G M + (positive energy) ^ trajecto- 
we have the following Hermiticity condition S*(X,k ) = [S/A*, k)}^ 1 which 
implies, as regard the scattering phases S of S/A, fc) = e 2 * <5 ( A > fc ) ; the relation 
S*(X,k ) = S(X*,k), so that they are real when A is also real. From what 
has been said above, S/A, k) is a meromorphic function in the variable A at 
least in the region 3ftA > 0, where /(A, ±k) are both analytic when k is real, 
so that, in the region 9ftA > 0, due to zeros of /(A, — k), poles of S/A, fc) 
may exist, which are universally known as Regge poles ; they are confined 
in the region SA > 0 and there exist independently of the potential U{r), 
that is, their existence has a pure kinematical nature. If one considers further 
restriction on the potential U (r), then the localization of the Regge poles may 
be more accurately identified. For instance, for purely imaginary values of r, 
if \U(ip)\ < Mjp 2 with M G M + , then the Regge poles are confined, in the 
half-plane 9ftA > 0 and TyA > 0, within the region delimited by the hyperbole 
branch SRA-SA < M/2 and by their asymptotes 3?A = 0 and 3A = 0. Instead, 
if \U(ip)\ < N/p with N G M + , then the Regge poles are confined into the 
region !RA < N/k, whereas, for fixed values of k G M + , the half-plane on 
the right hand side of the fine 9ftA = N/k, is a region in which S(X,k) is 
holomorphic. Moreover, for Yukawian potentials verifying \U{ip)\ < M/ p 2 , 
it is possible to prove that the real part of the Regge poles is upperly limited 
by 3ftA = M' so that the half-plane 7?A > M' is a holomorphic zone for 
S(X,k). We have also need to know the behavior of S(X,k ) as |A| —> oo in 
the region 9?A > 0. To this end, it is possible to prove that the following 
relation holds 


(127) lim S( A, k) = 1 

| A|—>oo 

at least into certain regions of the A plane. Such a relation may be also 
17 We have 1 RE' - E - it) = P/{E' - E) + in5(E' - E). 
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deduced by some estimates like the following ones. If S( A, k ) = e 2lS ^ x ' k \ then 

(128) \S(X, k) - 1| = 0(\- 1/2 e~ aX ) 

for potentials |f/(r)| < ce~ mr /r such that cosh a = l + m 2 /2/c 2 . For complex 
values of A, at least when the real part of the poles of S(X,k ) is upperly 
bounded, for |A| great enough and such that |arg A| < 7 t/ 2, we have 

(129) \S(X,k)-1\ <ce~ ayix , 

so that S(X,k) — 1 exponentially tends to zero as A —* oo in the region 
9ftA >0. It is desirable that many results of potential theory, when it is 
possible, may be extended to the relativistic case. In doing so, with respect 
to natural units, it is better to replace the fundamental kinematical variables, 
the relative wave vector k and the scattering angle 6, with the following two 
other variables, namely s = E = k 2 and the new parameter 

(130) t = q 2 = 2k 2 (l — cos#), 

which respectively represent the relative energy of the incident particle beam 
and the square of the transfer momentum during the scattering process. So, 
the scattering amplitude f(s,t), expressed in these two new variables, has a 
physical meaning only when s G M + and t G M _ fl ] — 4s, 0[. But, from these 
values, it is possible to analytically prolong such a function, as well as its 
asymptotic behavior, to all complex values of s and t, whose analytical prop¬ 
erties are often obtainable by means of certain dispersion relations involving 
such a function, like (124), which is a dispersion relation at a fixed 6 (hence 
at a fixed t), or the so-called double representation (in s and t) of S. Mandel¬ 
stam (see (Chew 1962)), which allows to determine the scattering amplitude 
without make use of any wave equation. The asymptotic behavior of f(s,t) 
as t —» oo has no physical meaning in the non-relativistic context but only in 
the relativistic one where the particle/antiparticle symmetry, which entails 
invariance with respect to the exchange between the two variable s and t, 
provides physical meaning to such a behavior. Indeed, at relativistic level, 
the so-called substitution law holds, according to which, together a diffusion 
process of the type 

(131) a + b —> a' + b', 

it is needed to consider the following cross-diffusion process 

(132) a + a' —» b + b' 
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where a denotes the antiparticle of a, and so forth. In a few words, the 
substitution law states that, in a diffusion process, an incoming particle (like 
a') is equivalent to its outgoing antiparticle (like a). The quantity t = q 2 , is 
the transfer momentum into the direct channel (131) as well as the energy 
of the cross channel (132). t is a limited quantity for the former, being 
0 < t < Ak 2 , whilst isn’t for the cross channel for which t may take any 
value in M + , so that it is possible to consider the asymptotic behavior of the 
scattering amplitude 


(133) 


f(k, 0) = f\y/E , arccos (l - — )j 


F(E,t) 


as t —> oo. If we work at a fixed value of energy, that is to say k e M + is 
arbitrarily fixed, then, by (130), it follows that t —> oo is equivalent to cosd —y 
—oo, so that it is needed to search for a representation of the scattering 
amplitude f(k,6), which is valid for physical values of 9 and that may be 
analytically prolonged in a open region of the plane of cos 9 in which it is 
possible bringing cos 9 to — oo. 

In pursuing this, we start from the following partial wave expansion of 
the scattering amplitude 


(134) /(M) = J]( 2/ + lMA;)P,(cos0) 

1=0 


which defines f(k,9) as an analytic function of cos 9 in that region of the 
plane of the complex variable z = cos 9 in which the series at the right hand 
side of (134), converges. Such a convergence region depends on the Z-th term 
of the series as l —> oo, so that we should study the behavior of ai(k) and 
Pi(cos9) as l —» oo. To this end, we have the following asymptotic relations 


(135) 



where a: is a positive real constant as above defined by cosh ck = 1 + m 2 /2k 2 , 
so that the Z-th term of the series (134) exponentially goes to zero as l —> oo 
in the region \ ( o9\ < a which, therefore, identifies the convergence region of 
the series (134) in the 9 plane. To determine the strip |9$| < a in the plane 
of z = cos 9, let us set x = y = ^sz, 9\ = 3 $9 and 9 2 = ^s9 , whence 


(136) z = x + iy = cos 9 = cos (9i + i9 2 ) = cos 9i cosh 9 2 — i sin 6\ sinh 9 2 


from which it follows that 


(137) x = cos 9\ cosh 9 2 , y = — sin 9\ sinh 9 2 , 
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so the boundary of the convergence region of (134) is the boundary of the 
strip |S0| < a, that is to say 9 2 = ^s9 = ±a, and inserting these values of 9 2 
into (137), hence eliminating 9 1 , we get the equation of the boundary of the 
convergence region of (134) in the plane of z = cos 0, that is the ellipse 

x^ ip 

(138) —^- + - 72 - = !. 

cosh a smh a 

said to be the Lehmann ellipse, (see (Lehmann 1958)) having major radius 
cosh a = 1 + m 2 /27c 2 , minor radius sinha = \J cosh 2 a — 1 and foci in the 
points cos 9 = ±1, the extremes of the physical region related to (134), which 
identifies the largest ellipse in the cos 6 plane where the scattering amplitude 
is analytic. It is possible to prove that, through (137), the strip |7y0| < a is 
the interior of the region identified by (138), so that the convergence region 
of (134), in the plane of cos 0, is the interior of the Lehmann ellipse. In 
the plane t = 2k 2 {l — cos 0), the Lehmann ellipse remains an ellipse, with 
foci in t — 0 and t = Ak 2 , the point t = —m 2 is the extreme of the ma¬ 
jor radius corresponding to cos 0 = 1 + m 2 /2k 2 , its two semi-axes depends 
on the energy through the above expressions of cosh a and sinh a, and the 
more the energy grows up, the smaller and the more restricted such an el¬ 
lipse becomes, at last reducing to the physical region given by the segment 
[—1,1] C JLs of the plane of z = cos 9 as k —> oo. Therefore, (134) is a repre¬ 
sentation of the scattering amplitude which converges within the Lehmann 
ellipse where it turns out to be an analytic function of cos 6. Nevertheless, 
this representation is unable to study the behavior of the scattering ampli¬ 
tude as t —> oo or equivalently as cos 0 —oo that, however, may be suitably 
re-written to accomplish this purpose transforming the series into an integral 
over the complex plane of the above considered variable A = l + 1 / 2 , in such 
a manner to have an analytic continuation of f(k, 6) over an open domain of 
the plane of cos 0 ; such a particular integral transformation is usually called 
Watson-Sommerfeld transformation. To introduce this, we need to consider 
the Legendre function Pi(z) = P\-i/ 2 (z) for arbitrary complex values of A, 
which is an even entire function of A 2 . For z = cos 9, we have the following 
estimate 

(139) |P A _ 1 / 2 (cos 0)| < c| sin0|- 1 / 2 |A| _ 1 / 2 e |SM,||9f ® l_9f ®' SiA . 

Now, we consider the function 

(140) G(X) = A <S(A,/c) ~ 1 P a _ 1 / 2 (- CO s0) 

COS 7T A 

which is a meromorphic function of A in 9ftA > o, with poles in the zeros of 
cos 7 tA (in correspondence to the physical values of A = l + 1/2) plus further 
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possible poles in 71A > 0 given by S( A, k ). The residue of G(A) in a generic 
pole in the variable A = / + 1/2, with / having physical values, is 

(141) Res{G(A)} A=m/2 = ~(21 + l)[S t {k) - 1]Pj(cos0). 

If one defines the integral 

(142) J G(X)dX 

where C is the path surrounding all and only all the poles of the real semi-axis 
9ftA > 0, then, from the residue theorem, we have 


(143) 


G(X)dX 


-2m ^ Res{G'(A)} A=z+ i /2 

1=0 


—2kf (k, 6) 


whence we have the following integral representation 
< 144 > '<*•«)= 

We would like to homotopically deform the path C in such a manner it 
coincides with the imaginary axis. First, we suppose do not exist poles of 
S(X,k ) in 5RA > 0. We notice that, as |A| —>■ oo along any direction of the 
half-plane 9ftA > 0, from (139) with cos 9 —» — cos 9 = cos(7r — 6), it follows 


(145) 


AP A _l/2(—COS0) 


<c(sind)- 1 / 2 |A| 1 / 2 e-l^ A l +oe - SRA , 


COS 7T A 

and, if ^(A, fc) — 1| < ce~ aXX , then, from (144), it follows that 
(146) |G(A)| < c |A| 1/2 e -l^- aA l + ( a0 - Q)5RA 


whence, in turn, it follows that G(A) goes to zero uniformly as |A| —> oo along 
any direction in JiA > 0, provided that |9?0| < a, which is a condition surely 
satisfied within Lehmann ellipse. Therefore, in these conditions (above all, in 
absence of poles of S( A, k)), the deformation of C to the imaginary axis SA is 
allowed, so obtaining the following integral representation of the scattering 
amplitude 


/(M) = 


(147) 


1 

2k 

1 

2k 


—too 

rioo 


G(X)dX = 
S{X,k ) - 

COS7tA 


-P\- 1 / 2 (— cos 9)dX 
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which is said to be the Watson-Sommerfeld transform (briefly WS transform ) 
of S(X,k ); it expresses the scattering amplitude f(k,9 ) as a continuous su¬ 
perposition of conical functions, that is to say, Legendre functions of the 
type Pa- 1 / 2 ( 0 ;) with t 6 I, Following (De Alfaro & Regge 1965, Chapter 
13, Section 13.2), the formula (147) looks like an infinite level Breit-Wigner 
formula. In (147), the variables A and cos 9 are conjugate of each other by 
means of the so-called Melher-Fock inversion formula through which it is 
possible to consider S(X,k) in terms of f(k,9 ); therefore, the properties of 
f(k,9 ) with respect to the variable cos 9, may be deduced by the properties 
of S(X,k) with respect to the variable A, and vice versa. Of course, further 
questions of convergence and related areas of convergence of the integral of 
(147) occur: in any case, what is important is to notice that there exists a 
common convergence region in which both representations (134) and (147) 
are the analytic continuation of each other; in particular, the representation 

(147) comprises the asymptotic region cos# — y —00, so that the Watson- 
Sommerfeld transform is suitable for studying the behavior of the scattering 
amplitude as cos 9 —> —00, that is, as t —» 00, whose properties in function 
of t — 2k 2 {l — cos 9) are obtainable from those in cos 9. 

It is interesting to consider the Watson-Sommerfeld transform of S( A, K) 
in presence of Regge poles, that is to say, in the case in which S(X,k ) has 
a certain number of poles in the quadrant I? A > 0, A A > 0. If, for each 
fixed value of k, S( A, k) has a pole of the type A = a n (k) + 1/2 with residue 
a n = cr n (k), then the residue of the function G(A), given by (140), computed 
in this pole, will be 

(148) {Res G(A)}A= a „+i /2 = — (a+ + 1/2) —— ^—^—P a (— cos#). 

Sill 7TQ: n 

Therefore, in applying the Watson-Sommerfeld transform along the deforma¬ 
tion of the original path C to the imaginary axis, we must further consider a 
closed path which surrounds every single pole of such a type when we apply 
(143) in pursing this, each of which provides the following contribution 


(149) 


Pn(k) 

sin 7T ot n {k) 


Pa n (k)(- COS 9) 


where 

? 7r I 

(150) /3 n (k) = — [a n (k) + -]cr n (k), 

18 See (Mehler 1881) and (Fock 1943), as regard the original papers. 
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so that, in presence of Regge poles, the Watson-Sommerfeld transform as¬ 
sumes the following expression 

9)d\+ 


the summation being over all the existent Regge poles. This is the funda¬ 
mental formula to study the asymptotic behavior of the scattering amplitude 
f(k, 9) = F(E, t ) as cos 9 —» — oo, that is, as t —» oo, which we now will carry 
out only in the case of a unique Regge pole A = a(k) + 1/2, in such a manner 
that (151) reduces to 


/(M) = 


(151) 


2k 


X S(A,*).-!„ , 

A - 7 - P\—l/2\ COS 


Pn(k) 


sin 7 ra n (k) 


COS7tA 


P«„(k)(~cos6) 


(152) 




no ° S(X, k) —1 

A-t—Pa- 1/2 (— cos 0)dA+ 


0(AO 


sin 7ra(A;) 


COS7tA 


P«(fc)(- cos 0) 


with k = a/P and — cos 9 = (t/2k 2 ) — 1 = (t/2E) — 1. Now, from the 
asymptotic behavior of Legendre functions, it is possible to deduce that the 
integral at the right hand side of (152) goes to zero as t —> oo, which is the 
main reason for which it is often called the background integral ; furthermore, 
it is also possible to prove that Regge pole behaves like t a ^ k \ so that, in the 
case of a unique Regge pole with Ra > —1/2, this dominates the background 
integral, so obtaining the following asymptotic expansion of the scattering 
amplitude given by (133) 

(153) F(E,t) ~ c(P, Res G(a(k)))t a{k) (Regge’s theorem) 


where c(E, Res G(a(k))) is a constant which depends on the energy and on 
the residue of the Regge pole. The asymptotic formula (153) still holds for a 
finite number of Regge poles in which a will be replaced by that pole a n with 
highest real part, whereas, in general, it does not hold for an infinite number 
of poles. Regge’s theorem is a consequence of Watson-Sommerfeld repre¬ 
sentation of scattering amplitude. Attempts to extend toward relativistic 
regimes - or, however, to high energy contexts - the relation (153) have been 
made: in this latter case, the (153), thanks as well to the substitution law, 
shows that the high-energy behavior of the scattering amplitude is mainly 
ruled by the related Regge poles involved in the given collision process, be¬ 
sides certain Regge cuts of the quadrant 9ftA > 0, QA > 0 whose existence 
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seems to be suggested by various physical models and that contribute as well 
to the asymptotic behavior of the scattering amplitude. We finally will try 
to shed light on the possible physical interpretation of the Regge poles. 

A single Regge pole provides the following contribution to the scattering 
amplitude, 

(154) — P a {— cos 9) = f R (k,0) 

Sill Tra 

where the position of the pole a(k) + 1/2 depends on the energy E = k 2 , 
while [3 is a residue which, in turn, depends on the energy as well. From (65) 
and (154), as well as taking into account certain properties of the Legendre 
functions, it follows that the contribution of this pole to the Z-tli partial 
amplitude is 

(155) af(k) = ^ f f R (k,0)P l (cos0)dcosd = n 

with a = a(E),f3 = /3(E), from which it turns out that this contribution 
of the pole to the Z-th partial amplitude is very relevant when a ~ l: for 
instance, if, whilst E is varying, a(E) takes, for E = E 0 , a value, say «o, much 
near to an integer one, say /, then the Z-th partial wave shows a resonance in 
correspondence at E — Eq. This last argument may be also supported by the 
fact that (155) is closely related to the Breit-Wigner formula (20) because 
it is possible to prove that it reduces to a formula of Breit-Wigner type. In 
fact, we suppose that S( A, k ) has a Regge pole for A = a(E) + 1/2 and that, 
for E = E 0 , let a(E 0 ) = a 0 be much near to an integer value l, that is to say 
a(E 0 ) = a 0 = / + ? 7 o + fcr 0 wher€@ |? 7 0 | -C 1, cr 0 -C 1. But, if S/A,/:;), for A; = k 0 
(that is to say, for E = E 0 ), has a pole for A = a(E 0 ) + 1/2, then the Jost 
function /(A, — k) is zero in these values, that is f(a(E 0 ) + 1/2, —Zc 0 ) = 0. 
In a neighborhood of E = E 0 , the position of the pole in the plane A will 
be given by A = A (E) = a(E ) + 1/2, where the function a(E) is defined 
as a solution to the equation f(a(E) + 1/2,—A;) = 0. But, since f(X,—k) 
is locally analytic in the variables A and k, with df /d\ = df /da ^ 0, it 
follows that a(E) will be, in turn, locally analytic in the given neighborhood 
of E = Eq, so that we may write a(E) l + tjo + zero + ct' 0 (E — Eq) where 
ccq = (da/dE) E=Eo = r/Q + ia' 0 and a' 0 ,rj' 0 real numbers. Therefore, in the 
neighborhood of E = E 0 , due to \r]' 0 \ <C 1, <r 0 -C 1, we have (Z + a + 1) ~ 2Z + 1 
and l — a ~ — ^[E — E 0 + ( i] 0 + zcr 0 )/ao] = —a' 0 (E — E 0 + A E + ?T/2) where 

19 The Regge pole a(Eo) = ao = (l + rjo) + iao , has positive imaginary part uq because 
every Regge pole lies in the region +A > 0. 
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A E = (/7o/7o+ <J o°‘o)/l a ol 2 an d r = 2(a 0 r]' 0 —r] 0 a , 0 )/\a' 0 \ 2 , so that (155) assumes 
the following expression 


(156) 


a ?{k) « 


1 _ P/a'p _ 

ir(21 + 1) E - E 0 + AE + *r/2‘ 


Under the hypotheses |r/ 0 1 <C 1 and <j 0 1, we also have AE <C E 0 and 

T <C E (h whence, in these conditions, the l -th partial wave basically will 
coincide with af, so that the contribution of the l -th wave, given by 


(157) 


oi(k) = 4tt( 2/ + l)|a ; (l-)| 2 , 


to the total scattering cross section, will be 


(158) 


„ 4\(3/aX/m + l)n] 

l{ J ~ (U - U 0 + AU) 2 + T 2 /4 


which is a formula of Breit-Wigner type, showing that, under the above 
hypotheses, we have a resonance, in the wave /, for an energy E = E 0 — AE, 
with width T. The number AE represents the small variation between the 
real energy E 0 , with which the Regge pole goes through nearest possible to 
the physical value l + 1/2 in the plane of A, and the real part Eq — AE of 
the energy with which the resonance appears in the plane of E, for l having 
a physical value. 

Now, as, in the E plane, the imaginary part - proportional to the width 
T - of the pole associated with a resonance, has a well-determined physical 
meaning related to the mean lifetime r = 1/T, likewise it is possible to give, 
in the A plane or plane of the angular momentum, a direct physical meaning 
also to cr 0 , i.e. the imaginary part of the pole associated with the same reso¬ 
nance, as follows. If one imagines, from a classical standpoint, a resonance as 
a (semi-bound) metastable system in which the related constituents rotate 
around each other for a time r, then <Jo is related with the angular lifetime 
of such a resonance. Indeed, if A = A (E) is a value of A for which a Regge 
pole appears, then we have 


(159) 


f(\(E),-k) = 0 


and the regular solution <j)(\(E), k(E),r) undergoes the asymptotic behavior 
given by (103), that, from (159), reduces to 


(160) 


2 ik 


f(\(E),k)e 


ik-f 
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If Afc > 0, we consider the following equations 


(161) 


where 


(162) 


0 " + 
0 " + 



A 2 (E) - 1/4 
r 2 

A 2 (E) - 1/4 


U(r))(/> 

U{r))i> 


0 


-0 + 


dX 2 (E) 0 
dE dE 


d(j)/dE. From them, it follows that 



so, when E and A are quasi-real number^, then also 0 2 will be a quasi- 
real and non-negative function, hence also d\ 2 /dE will be a quasi-real and 
non-negative function as well, to be precise equal to R 2 if R is the classi¬ 
cal radius of the orbit of the system, so that, if a' = da/dE = dX/dE = 
(l/2X)(dX 2 /dE), then we have a 1 ~ R 2 / 2A. Furthermore, Co is linked to the 
width T of the resonance by the approximate relation T/2 « do/ag, so that 
we have <7o ~ R 2 T /4A that, in natural units, reduces to cr 0 ~ RT/2v because 
the angular momentum approximately is A ~ (1/2 )vR. Therefore, if At is 
the revolution period, then we have At = 2 itR/v, and being r = 1/T, we have 
too do ~ (l/47r)(At/r), so that <Jq/2 represents the angular mean lifetime A 6 
of the system, that is 1/2(T 0 ~ 2n{r/At) = A6, by which it follows that the 
lower the imaginary part of the Regge pole is, the greater A 6 will be. Accord¬ 
ingly, we may conclude that a Regge pole, for A = a(E) +1/2 = cei(E) + 1/2, 
describes, at varying the energy E, a trajectory in the plane A; each time that 
a Regge trajectory on = ai(E) goes through very near to an integer number 
/, then a resonance will appear. Thus, the function cti(E) itself may give 
rise to many resonances, which may also be assembled into families, each of 
which associated with a single representing trajectory a = a n (E), this last 
idea having given rise to to many other applications and phenomenological 
speculations in the held of high and low energy physics. Therefore, to sum 
up following (Eden et al. 1966, Chapter 3) and (Mackey 1978, Section 21), 
bound states of Schrodinger equation for a spherically symmetric potential 
fall into families characterized by increasing angular momentum and decreas¬ 
ing binding energy, such a family appearing as a sequence of poles occurring 
in the successive partial wave amplitudes a/(s) l G No, at increasing values 
of s. As already mentioned above, following a previous analytic continuation 

20 See, for instance (Kontolatou 1993) and (Goodman & Hawkins 2014, Definition 2.2). 
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technique employed by S. Mandelstam in 1958 in the case of f(k , 9), the the¬ 
ory of complex angular momentum was introduced by Tullio Regge, in 1959, 
thanks to which it was possible to show that one could extend Si(k) in a nat¬ 
ural way as to be defined for all positive real l, and that therefore l —y Si(k ) 
could be regarded as a boundary value of an analytic function of two complex 
variables, the famous Regge poles being the poles of / —y Si(k ), with respect 
to l and for fixed real values of k, located in the half-plane dtl > —1/2 where 
Si(k) is meromorphic with respect to l. To be precise, the partial scattering 
amplitude ai(k) undergoes an analytic interpolation from integer valued to 
complex values that, for a certain class of Yukawian potentials, is uniquely 
determined by the asymptotic behavior of AZ —y oo, thanks to Carlson’s the¬ 
orem. At the same time, with this analytic continuation of Si(k), considering 
s — k 2 with its negative and positive determinations, we may also analyti¬ 
cally extend f(s,t) by means of a sum of partial scattering amplitudes ai(k) 
via the Watson-Sommerfeld transformation as given by (151) and the relation 
(64) which links together scattering phases (i.e., Si(k)) with partial scattering 
amplitudes (i.e., ai(k)). Therefore, the precise connection between asymp¬ 
totic behavior and individual angular momentum singularities was usually 
achieved by Regge’s method via Watson-Sommerfeld transform (see (Chew 
1966, Chapter 8, Section 8-2.)). Varying k G R, we obtain a family of curves 
in the l plane, say Si(k n ), hence a family of poles, say ai(k n ), which therefore 
change with k, so obtaining a discrete family of functions (Regge trajecto¬ 
ries) k —> o>i{k ) at varying of k, where cq(fc) is a pole of / —» Si(k) for each k 
fixed, which identifies its relative position on the given Regge trajectory. For 
s G M - , ai(k) is a real increasing function with respect to k] then, whenever 
oti{k) passes through a positive integral value Iq, then the corresponding value 
of k is a pole of Si Q (k ) corresponding to a bound state with negative energy; 
for k G M + , cq(s) is instead a complex function, with positive imaginary part, 
so that when k passes from negative values to positive ones, the Regge tra¬ 
jectory leaves the real axis and enters into the upper half-plane 3Z > 0 of the 
complex l plane. Therefore, the bound states of the dynamical system lie on 
Regge pole trajectories, several times on each trajectory, so that the bound 
states may be grouped into families, all members of the same family lying on 
a single Regge trajectory. Again, following (Eden et al. 1966, Chapter 3), 
bound states of Schrodinger equation for a spherically symmetric potential 
fall into families characterized by increasing angular momentum and decreas¬ 
ing binding energy, such a family appearing as a sequence of poles occurring 

21 Because such a function is initially defined, in the variable l , only for an infinite discrete 
set of values, that is No, so that it is more correct to speak of an analytic interpolation 
rather than an analytic continuation. 
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in the successive partial wave amplitudes a*(s) l G No, at increasing values of 
s. The Regge theory pictures this sequence as due to the presence of a single 
pole whose position varies continuously with l and which is relevant to the 
physics of bound states only when l takes non-negative integral values. This 
idea has given meaning by the construction of an interpolating amplitude 
az(s), defined for non-integral, and indeed complex, values of /, which coin¬ 
cides with the physical amplitudes ai(k) when l G No- This function oj(s) is 
an analytic function of its arguments except for certain singularises. Among 
these singularities will be a pole (Regge pole) corresponding to each of the 
bound state families, the location, or trajectory, of such a pole being given 
by the equation l = a(s), the bound state energies corresponding to values 
of s which make / take the values 0,1,2,.... One can deal similarly with 
unstable bound states and resonances: indeed, given any pole k 0 of oii 0 {k) in 
the complex k plane, we may obtain a nearby pole in the complex k plane 
for each value of l near / 0 . Thus, we obtain a two-dimensional submanifold of 
the direct product of the complex / plane with the complex k plane, having 
the property that each pair (/, k) in this submanifold is a pair such that k is a 
pole of ai{k). These submanifolds will be finite or countable in number, and 
we may group together resonances when they lie on the same one, speaking 
of Regge recurrences of a given bound state or resonance. To be precise, if 
there exist positive values of k in which ai(k) has an integer real part and 
a very small imaginary part, say a(k) = L + ia'(k ) with 0 < a'(k) <C 1, 
then a resonance with angular momentum L and energy k, does appear. The 
greater the value of k, the more the Regge trajectories depart out from the 
real axis, towards the plane dtl > 0, which, for attractive potentials, form 
clusters or families of bound states and resonances. 

Therefore, to summarize, the connection between singularities of the scat- Duality 
tering amplitudes, considered as a complex function analytically prolonged 
with respect to the energy or the angular momentum, and the presence of sta¬ 
ble (bound states) or unstable (resonances) particles, played a fundamental 
role in those theoretical approaches of elementary particle physics in which a 
central role is given by the scattering amplitudes and their analytical proper¬ 
ties; moreover, the phenomenology of the 1950s high-energy physics plainly 
recalled the attention on the importance of singularities of S matrix the¬ 
ory with their physical meaning, pointing out on the notable contribution of 
bound states and resonances till to suppose that non-resonant contributions 
could be described by means of suitable averages of resonances. This last 
idea is the conceptual key of the so-called duality, which, in turn, was at the 
basis of the string theory, whose early origins may be retraced in the late 
1960s works of Gabriele Veneziano, and whose main study’s object is the so- 
called relativistic string, a concept generalizing that of particle associated to 
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a given quantum field, understanding this as represented by the vibrational 
modes of certain geometrical structures, having a well-defined dimension, so 
providing a new theoretical framework for the strong interactions whose phe¬ 
nomenology of the time shown a prevalence of either diffractive phenomena 
and resonances which, therefore, entailed the consideration of a mathemat¬ 
ical structure generating infinite states, like a vibrating string. Following 
(Marchesini et ah 1976), at the beginnings, attempts to apply quantum held 
theory for treating strong interactions were tempted, going beyond the initial 
and limited Yukawian work on the mesonic held. The EM analysis of the 
structure of hadrons has moreover identihed a complex structure formed by 
more elementary components, called partons, to which it will be allowed to 
apply possibly the ordinary quantum theory of helds. Therefore, from this 
standpoint, hadrons should be seen as bound states of partons, while the 
collisions between hadrons should be considered as collisions between bound 
states. Nevertheless, the resulting theoretical framework presented difficult 
formal problems like the one regarding the confinement of partons within 
an hadron, so that a necessary way to follow for overcoming this type of 
problems consisted in finding non-conventional held theories or even com¬ 
pletely new formal sights, like the one which revolves just around S matrix 
theory, of which one of the most promising formal model was that based on 
the above notion of duality as meant by Veneziano. If we consider a general 
exchange procesJ^l of the type a + b —» a' + b', and assume to be valid the 
hypothesis of duality, then the related scattering amplitude may be given 
by the sum of two main contributions: a low-energy contribution given by 
the Breit-Wigner formula (20), say af w , related to a decay of resonances of 
the type a + b —> (R ) — y a' + b', and an high-energy contribution given by 
the Regge behavior according to (153), say af ~ ( k 2 ) a , so that we should 
have ai ~ af u + af. The phenomenological analysis of such a question has 
entailed that, at least in the average approximation, the related imaginary 
parts should be equal amongst them, that is 7 sai = ^saf = 7 saf w , so that, 
roughly, the Regge behavior and the formation of resonances basically are 
dual descriptions of the same phenomenon. Following (Martin & Collins 
1984, Chapter 7), one of the main characteristic properties of hadron inter¬ 
actions is that the scattering amplitude may be given as a sum of resonance 
contributions of the type a; ~ ^ res af" with af " given by the Breit-Wigner 
formula. In the limit of small resonance width T, we may also write this last 
sum as a sum of resonance contributions of mass rrij along the direct channel 

22 It is a process in which, for instance, a condition of the following type holds, i.e., the 
quantum numbers of a and a' are different between them; this is the case which happens, 
for example, when the exchange of quarks and gluons between the hadrons is involved (see 
(Collins & Martin 1984, Chapter 7)). 
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as follows 


(163) ai(s,t) = ^ 2 Cj(i) 

where s — E 2 — k and t is defined as above and proportional to the square 
of the four-momentum. But, in the limit T —» oo, from what has been said 
above, we have that the scattering amplitude may be also written as a sum 
of Regge contributions as follows ai ~ )T) res so that, along the crossed 
channel, we have too 

(164) ai(s , t) = V 1 2 Ci(s), 

t — mf 

i L 


so that duality implies that the scattering amplitude may be equivalently 
written both as a sum of resonances related to the direct channel (equation 

(163) ) and as a sum of resonances related to the crossed channel (equation 

(164) ), from which it follows another of the most characteristic properties of 
hadron interactions, that is to say the crossing symmetry in the variables s 
and t (duality), so that each of these two formulas expresses a single equiva¬ 
lent contribution to the scattering amplitude, this last property having been 
elegantly expressed by the mathematical properties of following formula for 
the scattering amplitude 


(165) 


ai(s,t) 


r(l -q(g))T(l -a(t)) 

T(1 — a(s) — a(t) 


proposed by G. Veneziano in 1968, where a(x) = a 0 + a'x is a linear 
Regge’s trajectory, T is the Euler gamma function with poles in No such that 
r(z) « 2 d. The symmetries of (165) take into account the duality principle 
of above; furthermore, it has poles of resonance type for s = m?(i — ao)/a', 
with highest spin cq(?7?|), that, for great values of s, gives rise to the Regge’s 
asymptotic behavior ai(s,t ) ~ s a ^\ The interest of (163) and (164) is that, 
the duality concept underlying them, may be extended to the case of scatter¬ 
ing amplitudes related to an arbitrary number of particles when one is able 
to identify the right variables generalizing s and t, so that such a duality 
hypothesis (with T —* 0) determines the form of every hadronic scattering 
amplitude, so casting the bases of a theoretical framework of hadronic inter¬ 
actions. Nevertheless, not always the basilar condition T —» 0 has a reality 
meaning because, in general, T grows together with the resonance mass. 

To this point, it is maybe better to follow, almost verbatim, (De Alfaro 
et al. 1973, Chapter 9), where a new approach to strong interactions is 
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presented just making use of the notion of duality. In (De Alfaro et al. 1973, 
Chapter 9, Section 1), the authors point out that in the recent theoretical 
development classified as ’duality’ the substitution rule, or ’crossing’, plays 
a fundamental role. Roughly speaking, in a typical process 

(166) A + B -> C + D, 

the same scattering amplitude should describe the process (166), together 
with the ’crossed’ processes 

(167) A + C -»■ B + D, A + D^ B + C. 


In the case of amplitudes with a larger number of external lines, then the 
number of processes which should be simultaneously described increases very 
rapidly with the number of external particles. It should be emphasized that 
the very strong predictive power of crossing is present only if one makes fur¬ 
ther assumptions amongst which analiticity and ’good’ asymptotic behavior. 
Indeed, since, for example, the three physical regions - i.e. the allowed ranges 
for the kinematical variables - of the processes (166-67) are separated and 
disconnected, we should carry out three independent descriptions of these 
three processes. Instead, the addition of analyticity assumptions allows, for 
instance, a continuation from one physical region to another, so enabling 
us to obtain strong constraints among these three different descriptions. In 
fact, it is very probable that an approximate amplitude which is perfectly 
reasonable for reaction (166) proves to be roughly inadequate for the ’crossed’ 
channels given by the reactions (167). So, the addition of the requirement 
of ’good’ asymptotic properties, like those provided by Regge prescriptions, 
has the consequence that any description which is complete for one of the 
three channels, should also be complete for the other two. Let us consider, 
for example, an amplitude satisfying the dispersion relation 


(168) 


A(s,t) 


1 

7T 


7L4(s', t) 
s' — 


where it is assumed that the dispersion relation converges for a sufficiently 
spacelike value of t, and possible singularities in the third variable s are, for 
simplicity, neglected. The latter equation (168) shows that a knowledge of 
the discontinuity of A(s,t), related to those physical processes in which s is 
the energy, leads, in principle, to a complete knowledge of A(s, t) and also of 
the analytic properties in the variable t. For instance, the possible poles in t, 
related to possible particles exchanged in the t channel, should not be added 
in an ad hoc manner to the right-hand side of equation (168), but should 
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come from the divergence of the dispersion integral in the variable s'. In this 
regard, it is well-known that the Regge’s asymptotic law 

(169) 9h4(s,f) ~ (3(t)s a V 

s—> oo 


gives rise just to this type of behavior, so we see a direct relation between s- 
and t-channel descriptions of scattering. A more precise and very convenient 
form of expressing such interrelations is given by the so-called finite energy 
sum rules (FESR) of the form 


(170) 



s n ^sA(s,t)ds = (3(t) 


j^a{t)+n+l 

aft ) + n + 1 


which gives a direct relation between the s-channel description of scatter¬ 
ing (that can be inserted into the left-hand side of this equation) and the 
t-channel Regge poles appearing in the right-hand side. In particular, if 
we wish to evaluate the left-hand side of equation (170), inserting only the 
contribution of the resonance exchanged in the s channel, then the relation 
(170) leads to a very important duality property which, broadly speaking, 
can stated in the following form: 


«the t-channel Regge-pole description of the process corresponds to an 
appropriate average of the amplitude computed in terms of the s-channel res- 
onances». 


The duality standpoint has very important implications in the phenomeno¬ 
logical analysis of data: indeed, it tells us that ’interference’ models based on 
the coherent sum of Regge plus resonance contributions should be discarded. 
From the theoretical point of view, the duality relation between s-channel res¬ 
onances and f-channel Regge poles, strongly suggests the possibility of setting 
up a duality bootstrap. Indeed, the f-channel Regge pole interpolates a family 
of f-channel resonances, so that, if we are working in a completely crossing 
symmetric problem (such as, for example, nu —> rnr), then we arrive at a 
large number of self-consistency constraints for resonances exchanged in the 
process. A first semi-phenomenological treatment shows that adjustments 
can be made with straight line trajectories and when an infinite number of 
parallel daughter trajectories, namely a n (s ) = ct(s) — n, are also present. 
The most direct and explicit way of carrying out the duality-resonance boot¬ 
strap leads directly to the Veneziano approach, which was inspired by the 
use of the phenomenological approach based upon FESR and duality boot¬ 
strap programme applied to the analysis of the process ttu —> 7T7t, so giving 
rise to first dual-resonance amplitude satisfying crossing and Regge behavior. 


47 



Veneziano approach is therefore based on the idea that a complete, although 
approximate, description of scattering can be obtained in terms of s-channel 
resonances only. If we insert in the dispersion relation (168) the resonance 
contributions coming from the particles contained in the trajectory a(s) and 
in its daughters, we obtain 


(171) 


A(s,t) 


£ 


Cn (t) 
a(s) — n 


where c n (t) is a polynomial of order n in t. The crucial point is that the 
representation (171) is complete and should contain (in an implicit manner) 
the poles in the crossed variable t. Those poles should arise from the possible 
divergence of the sum on the right-hand side of equation (171), which should 
thus contain an infinite number of terms. We see that one can hope to achieve 
a purely resonance bootstrap only in the case of infinitely-rising trajectories. 
A precise mathematical solution to duality constraints, namely 


(172) 


A(s,t) 


oo 


£ 

n=0 


a(s) — n 


v - "' c n (s) 


has been given by the now famous beta function amplitude of Veneziano, i.e. 


(173) 


A (s,t) 


r(-q( S ))r(-q(t)) 

T(-a(s) — a(t) 


It is easy to show that such an amplitude satisfies the bootstrap require¬ 
ment given by equation (172) and, moreover, is the Regge-behaved both in 
the s- and £-channels. We can furthermore see that the Veneziano model is 
physically acceptable if all trajectories are straight lines (corresponding to 
infinitely narrow resonances) with universal slope. If one moreover compares 
this result with Regge asymptotism, it is therefore found that elementary 
crossing allows only one kind of trajectory, the straight line, while devia¬ 
tions from the straight line will come together with more complicated types 
of singularities like Regge cuts. Starting from Veneziano amplitude, it has 
therefore been possible to develop a general dual-resonance model, which en¬ 
abled to construct, in a consistent and systematic way, dual amplitudes for 
any inelastic process. For example, a formalism based on harmonic oscillator 
operators plays a fundamental role in this theoretical approach. But, for a 
more general and deeper view of duality theory of strong interactions, we 
refer to (De Alfaro et al. 1973, Chapter 9). 
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2. On some historical aspects of 1959 Tullio Regge paper 

To sum up, following (Collins & Squires 1968, Introduction), although the 
generalization of angular momentum, say /, to a complex variable, and the 
subsequent transformation of the partial-wave series to an integral in the /- 
plane, was therefore introduced by Watson in 1918 and revised by Sommer- 
feld in 1949, it was not until the pioneering work of Regge in 1959, within S 
matrix theory framework, that the possible usefulness of the idea in particle 
physics was realized. Regge showed that, for a wide class of potentials, the 
only singularities of the non-relativistic scattering amplitude in the /-plane 
were poles which move with energy, i.e. poles at positions / = a(s). These 
are now called Regge poles, and a(s) is called a Regge trajectory. The method 
was used by Regge to close an important (if not physically very interesting) 
gap in previous proofs of the Mandelstam representation for potential scat¬ 
tering. However, the fact that the Regge trajectory correspond to physical 
particles (or resonances) when a(s) equals an integer for positive s, and also 
determine the high energy behavior of the cross-channel amplitudes, claimed 
the attention of elementary particle theorists, and since the early 1960s, the 
subject had a vigorous, controversial and exiting history. In (Frautschi 1963, 
Chapter X), it is also referred that Regge transformed the scattering am¬ 
plitude to a new representation involving complex / just by means of the 
so-called Watson-Sommerfeld transformation. This method, says Frautschi, 
had a history stretching over several decades, and was used to study rain¬ 
bows, propagation of radio waves around the Earth, and scattering from 
various potentials. Regge’s original contribution lay in understanding the 
special features of complex angular momenta for scattering from superposi¬ 
tions of Yukawa potentials (i.e., the type of potential deemed to be relevant 
for relativistic scattering) and calling these features to the attention of high- 
energy physicists. On the basis of what told by E. Guth, Frautschi says that 
a related transformation similar to that handled by Regge, was previously 
used by Poincare and Nicholson in the early 1910s in connection with the 
bending of electromagnetic waves by a sphere, hence explicitly introduced 
by Watson in 1918 and later resurrected by Sommerfeld in 1949. In this re¬ 
gard, Frautschi also quotes a previous related work of Nicholson dating back 
to 1907, and not cited by other references. In any event, already Poincare, 
Nicholson and Watson were aware of the existence of certain poles in the 
scattering amplitude of acoustic and radio wave propagation problems, then 
systematically treated by Sommerfeld who provided a first exhaustive the¬ 
oretical framework mathematically similar to the one then established by 
Regge in dealing with De Broglie’s waves, whose pioneering works of the late 
1950s opening, for the first time, the use of such formal techniques in elemen- 
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tary particle physics, clearly underscoring differences and possible analogies 
between the acoustic-electromagnetic case and the quantum one. 

In the monograph (De Alfaro & Regge 1965, Chapter 1, Section 1.2), the 
authors give a compendious historical outline of potential theory. In par¬ 
ticular, they point out what usefulness has their approach in reconstructing 
the scattering amplitude via Born approximation method, which uses tech¬ 
niques based on the introduction of a unique analytic interpolation of partial 
scattering amplitudes for complex values of the angular momentum, follow¬ 
ing 1959 Regge work. Then, in the next Chapter 2, the authors summarize 
the main mathematical tools used throughout the text, amongst which is a 
Poincare’s theorem already quoted in the previous section 1, that they say to 
play a key role in potential scattering. This theorem concerns with general 
linear second order homogeneous differential equations of the form 

{17 A) -^ +p ( x )J±± + q ( x )f( x ) = o 

which, following (Whittaker & Watson 1927, Chapter X, Section 10.1), can 
be written in the form 

(175) 
where 

(176) f(x) = (p{x) exp {~\J P(v)dy ) 
and 

Tl s / \ 1 dp(x) 1 2 

(177) J{x) = q{x)---^---p{x). 

For example, Schrodinger equation for partial waves may be written in the 
form (175). The above Poincare theorem (see (Poincare 1884a,b)) has to 
do with equations of the form (175) depending on a parameter 77 through a 
function J of the type J{x, r]), supposed to be an analytic entire function of 
77 . Take now a solution ijj{x) defined by a boundary condition independent 
on 77 in an ordinary point P{x = c). The theorem states that ^(cc, 77 ) for 
fixed x , as a function of 77 , is also an entire function. The condition that P 
should be an ordinary point, can be relaxed provided the boundary conditions 
are still 77 -independent. Hence, De Alfaro and Regge say that many of the 
theorems which will be established later in the text, are really generalizations 
of Poincare’s theorem, this fact witnessing the importance of such a result 


d 2 cp(x) 

dx 2 


+ J{x)ip{x) = 0 
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in the formal framework of potential scattering. For instance, we have used 
Poincare’s theorem to prove analyticity properties of a regular solution to 
(71). In (Regge 1959), it was considered, for the first time, the possibility 
of introducing the angular momentum as a complex variable, showing the 
convenience of this procedure in proving that Mandelstam representation 
for potential scattering holds. Regge starts with the consideration of the 
Schrodinger equation wrote in the form 

(178) r - (£zm - 1 + U(x))l> = 0 

having chosen dimensionless variables, by putting x = kr , where r is the dis¬ 
tance from the origin, and k the (fixed) wave number. Here A is a generalized 
complex orbital momentum, and when it assumes positive half-integer values 
(hereafter referred to as the physical values), then we shall write A = j +1/2. 
The equation (178) is even in A, and further restrictions are assumed for U (. x ) 
(see (72)). Then Regge considers a solution to (178), say F(\,r),x), whose 
parameter r/ may vary in such a manner to comprehend other solutions to 
(178) (see (Regge 1959, (1.3) and (1.4)); moreover, a slight generalized form 
of the above mentioned theorem of Poincare, states that F(X, r],x), for fixed 
r) and x, is an entire function of A. At the point 6. of his paper, Regge 
establishes some results in the field of dispersion relations. As is well-known, 
these relations are statements of analyticity of the scattering amplitude as 
a function of the energy and of the transmitted momentum. Although the 
energy is kept fixed, it is still possible to derive, for special classes of poten¬ 
tials, enough properties as to guarantee for the existence of such relations. 
As we have seen in the previous section, analyticity in cos 9 is known to 
subsist within the Lehmann ellipse (138), whose corresponding representa¬ 
tion is unable for studying the behavior of the scattering amplitude letting 
t —» oo or equivalently as cos 9 —» oo. But Regge cleverly and ably overcomes 
this formal problem transforming the series of the partial wave expansion of 
the scattering amplitude given by (134) (i.e., the expansion (6.3) of (Regge 
1959)), into an integral over the complex plane of the variable A = l + 1/2 
in such a manner to have an analytic interpolation of the scattering ampli¬ 
tude over an open domain of the cos 0-plane, so obtaining the integral (147) 
(i.e., the integral (6.4) of (Regge 1959) or the integral (9.19) of (De Alfaro 
& Regge 1965, Chapter 9, Section 9.3)), with an artifice that Regge says as 
due to G.N. Watson and used by A. Sommerfeld in some wave propagation 
problems, but without obviously giving other historical informations in this 
regard. In (De Alfaro & Regge 1965, Chapter 13, Section 13.2), the authors 
say that, strangely enough, formulas of the type (147) for general singular 
potential were first historically treated in connection with the propagation of 
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waves all around the Earth, where the potential U (r) was a truly hard-core 
potential, i.e. U (r) = oo for r < R, where R stood for the radius of the Earth 
supposed to have a spherical form besides to be perfectly conducting. The 
resulting pole expansion was seen to converge very rapidly for large k but 
there is no indication that this result may be extended to cases of significance 
in nuclear physics, as deemed by De Alfaro and Regge in the 1960s. 

Because of the importance of this last crucial formal passage with respect 
to the whole Regge’s work achieved in (Regge 1959) - whose primary aim 
was to prove the double dispersion relation in potential scattering - as well 
as for the notable role played by it in the next developments of high-energy 
physics, we now briefly devote a lot of time to deepen just these latter histor¬ 
ical aspects of such a 1959 Regge’s seminal paper. Following (Newton 1964, 
Chapter 1), the mathematical analysis of the results of scattering experiments 
has lately acquired a new tool based on the replacement of an infinite series 
of discrete terms by a contour integral of an analytic function introduced by 
Poincare and Nicholson, and first applied in its present form in electromag¬ 
netic theory by Watson in 1918 and, after its resurrection, in electromagnetic 
wave problems, by Sommerfeld. Also in (De Alfaro & Regge 1965, Chapter 
9, Section 9.3) some further historical clarifications are included. Indeed, the 
authors recall that the form of the scattering amplitude (147) is essentially 
derived by the work of Poincare in diffraction of electromagnetic waves and 
exposed in (Poincare 1910) which, nevertheless that, is not quoted in the 
original Regge’s paper of 1959. In (Poincare 1910), the author considers ap¬ 
plications of this his theorem quoted above to certain differential equations of 
the type (178) (like, for example, the equation (2) of Section 12 or the equa¬ 
tion (1) of Section 11), as well as Legendre polynomial expansions similar to 
(134) (like, for example, the series (2) of Section 6) and integral expression 
similar to (147) (like, for example, the expressions (3) or (13) of Section 13 
or the expression (3) of Section 15; see also Sections 17 and 19). As his 
usual style, Poincare does not quote any reference, except a mention to Max 
Abraham (1875-1922), at page 170, and to Arnold Sommerfeld (1868-1951), 
at page 171. Afterwards, De Alfaro and Regge quote too the works of John 
William Nicholson (1881-1955), namely (Nicholson 1910a,b; 1911), which are 
centered around mathematical problems inherent in diffraction phenomena 
of electromagnetic waves sent by an Hertzian oscillator round the surface of 
the Earth considered as a sphere of perfect conductivity, taking into account 
the previous work made by Poincare on this subject in which the electromag¬ 
netic forces involved are expressed as an integral of Fredholm’s type which 
will be then expanded into a series of zonal harmonics. In Nicholson’s works, 
the main formal aspect of his treatment of the related physical problem is 
just centered on the technique of replacing a series with an integral extended 
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to certain infinite regions (see, above all, (Nicholson 1910b) where integral 
expression of the magnetic forces are deduced from their series expansions) 
whose values are then determined by means of the calculus of residues. On the 
other hand, as early as 1900s, in many investigations of applied mathemat¬ 
ics, it is necessary to sum a series into each term of which a large parameter 
enters. For example, many problems of diffraction of waves by large obsta¬ 
cles depend for their ultimate solution upon the sum of a series of harmonic 
terms in which complicated derivates of Bessel functions of large argument 
are present. By expressing these transcendental functions in an appropriate 
manner, these summations may all be made to depend upon a single type, 
with which the paper (Nicholson 1907) proposes to deal, and in which series 
to be summed are, at a first approximation, replaced by a suitable integral 
built up according to a certain new method which will be further retaken by 
Nicholson himself in the next papers of 1910s. Then, De Alfaro and Regge 
refer that such a mathematical technique was extensively used by George 
Neville Watson (1886-1965) in (Watson 1918) (see also (Watson 1944) and 
references therein), where Poincare’s and Nicholson’s works are quoted, to¬ 
gether that of H.M. Macdonald (see (Macdonald 1914)), H.W. March, M. 
von Rybczynski and A.E.H. Love (see (Love 1915) and references therein). 
In any case, ever since the time, about 1902, when G. Marconi first succeeded 
in sending wireless signals across the Atlantic ocean, the question of explain¬ 
ing the mechanism of such transmission attracted the attention amongst 
mathematicians and physical-mathematicians, many of whom are quoted in 
(Love 1915), amongst whom are H.M. Macdonald, J. Zenneck, H.W. March, 
M. von Rybczynski, L.W. Austin, J.L. Hogan, and W.H. Eccles (see (Love 
1915) and references therein), and this ever since early 1900s. But, just 
following what is said in (Love 1915, No. 10), for the infinite series which 
represents the effect of curvature without resistance, methods of summation 
were devised by Macdonald (see (Macdonald 1911)), Poincare and Nicholson 
(see their above quoted works, in addition to (Poincare 1903)); and again 
another method was devised by Macdonald (see (Macdonald 1911)). All 
these methods depend upon a transformation of the series into a definite in¬ 
tegral, hence by an approximate evaluation of the latter integral. Poincare, 
in (Poincare 1910), did not press his method so far as to tabulate numerical 
results, but concluded that the expression for the electric force normal to 
the surface, at an angular distance 9 from the originating doublet, should 
contain a factor of the type e~ Ad . Nicholson, in (Nicholson 1910a,b; 1911), 
went further in the same direction, obtained a formula for the magnetic force 
containing such an exponential as a factor, and deduced definite numerical 
results. Macdonald, in (Macdonald 1911), also obtained definite numerical 
results which cannot be reconciled with those of Nicholson. The discrepancy 
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was discussed by Nicholson himself, who traced it to an alleged flaw in the 
analysis used by Macdonald in (Macdonald 1911); it was also discussed by 
Macdonald in (Macdonald 1914), who pointed out a difficulty in the analysis 
used by Poincare in (Poincare 1910) and by Nicholson in (Nicholson 1910a,b; 
1911). Fresh numerical results were then deduced by Macdonald in (Macdon¬ 
ald 1914) from a new method of summing the series, but they did not agree 
with those previously found by Nicholson, or with those previously found by 
Macdonald himself in 1911. Nevertheless, many of the authors and related 
works quoted in (Love 1915), which have to do with the above crucial passage 
from a series to an integral, are not yet mentioned in (De Alfaro & Regge 
1965). 

The 1918 paper of Watson, however, put into reciprocal comparison 
mainly the works of Nicholson and Poincare on the one hand, with the works 
of Macdonald on the other hand, so reaching to a major clarification either in 
physical-mathematics consequences and formal aspects of the original physi¬ 
cal problem under examination. Under advice of B. Van der Pol, who asked 
him to further analyze the problem with the chief purpose to clarify why 
there subsist the various discrepancies found by different authors in dealing 
with this physical problem, Watson first highlights, from a further compari¬ 
son between H.W. March and M. Von Rybczynski works with that of A.E.H. 
Love and with those of Poincare and Nicholson, what follows (see (Watson 
1918, No. 2)) 

«The essential advance in this paper is closely connected with the funda¬ 
mental error of March and Rybczynski which was pointed out by Love. In 
dealing with an oscillator on the positive half of the axis of harmonics, those 
writers express a Hertzian function by an integral of P s (cos 9), the integration 
being carried out with regard to the degree s of the Legendre function; such an 
integral has a line of singularities along the line 6 = n, and is regular along 
the line 9 = 0. The fact is that, when harmonics of non-integral degree are 
introduced, the appropriate function to use is not P s (cos9) but P s (— cos9); 
this fundamental point is somewhat obscured by the equation 

P n (— cos 9) = (-l) n P n (cos6>), 

which holds between the functions whose degrees are integers. The failure of 
convergence of an integral involving P s (cos9) along the line 9 = 0 (when an 
oscillator is placed on the positive half of the axis of harmonics) is strictly 
analogous to the failure of convergence of the series 1 + z + z 2 + all round 
the circle \z\ = 1 on account of the singularity of the function 1/(1 — z) at 
the point z = 1. A simple electrostatic example is afforded by the potential of 
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a unit charge at distance a from the origin. The potential near the origin is 


.. OO 

(*i) V ' = i E(-)” P »( '0=9) = 

a \aJ 

n=0 

n =o 

(* 2 ) = f- l (-)>.(cos9)— 

ha J \a / sms7r 

where the contour starts from +oo and returns to +oo a/ter encircling the 
points s = 0,1, 2,... which are poles of the integrand. On swinging round the 
contour so as to surround the other poles of the integrand, and evaluating 
the residues, we find the series for V in descending powers of r, valid when 
r > a». 


Watson, in discussing the swinging of the contour of integration of (* 2 ), 
makes reference to a paper of Ernest William Barnes (1874-1953) (see (Barnes 
1908)) in which the so-called Mellin-Barnes integral is defined, as well as he 
quotes a Laplace’s formula for P s (cos9) to prove the convergence of the inte¬ 
gral for |s| large (see (Watson 1944, Chapter 6, Section 6.5)). Therefore, it is 
more probable that Watson, in achieving the integral (* 2 ), considered what 
Barnes made in this respect in defining his integral which is, roughly, a con¬ 
tour integral of products involving gamma function and exponential factors 
(see also (Barnes 1910)). On the other hand, recent history of mathematics 
research tells us that it was Salvatore Pincherle (1853-1936), in his works 
on generalized hypergeometric functions, to have early applied the elements 
of the technique of contour integration methods to special functions (to be 
precise, to hypergeometric functions) with the introduction of first forms of 
integrals of Mellin-Barnes type. Indeed, following (Mainardi & Pagnini 2003, 
Section 1), 

«In Vol. 1, p. f9 of Higher Transcendental Functions of the Bateman 
Project, we read ”Of all integrals which contain gamma functions in their 
integrands the most important ones are the so-called Mellin-Barnes integrals. 
Such integrals were first introduced by S. Pincherle, in 1888 [see (Pincherle 
1888; 1965)]; their theory has been developed in 1910 by H. Mellin [see (Mellin 
1910)] (where there are references to earlier work) and they were used for 
a complete integration of the hypergeometric differential equation by E. W. 
Barnes [see (Barnes 1908)]. In the classical treatise on Bessel functions by 
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Watson [see (Watson 1944, Chapter 6, Section 6.5)] we read: ”By using in¬ 
tegrals of a type introduced by Pincherle and Mellin, Barnes has obtained 
representations of Bessel functions which render possible an easy proof of 
Rummer’s formula [...]. Here we point out that the 1888 paper (in Italian) of 
S. Pincherle on the Generalized Hypergeometric Functions led him to intro¬ 
duce the afterwards named Mellin-Barnes integral to represent the solution 
of a generalized hypergeometric differential equation investigated by Goursat 
in 1883. Pincherle’s priority was explicitly recognized by Mellin and Barnes 
themselves, as reported below. In 1907 Barnes [see (Barnes 1907)], wrote: 
’’The idea of employing contour integrals involving gamma functions of the 
variable in the subject of integration appears to be due to Pincherle, whose 
suggestive paper was the starting point of the investigations of Mellin (1895) 
though the type of contour and its use can be traced back to Riemann”. In 
1910 Mellin [see (Mellin 1910)], devoted a section (§ 10: Proof of Theorems 
of Pincherle) to revisit the original work of Pincherle; in particular, he wrote 
’’Before we are going to prove this theorem, which is a special case of a more 
general theorem of Mr. Pincherle, we want to describe more closely the lines 
L over which the integration preferably is to be carried out”». 

Finally, Sonnnerfeld, in (Sommerfeld 1949, Appendix II to Chapter V, and 
Appendix to Chapter VI), retakes the above mentioned work of Watson in 
computing the solution u to the wave equation A u + k 2 u = 0, of which a first 



form is obtained as a series expansion of the 




whose convergence, however, is very poor, and limited to the domain given 
by a < r < oo and 0 < 6 < n. This series expansion has been stud¬ 
ied either by P. Debye in 1908-09 as well as by P. Frank and R. von Mises 
in the second volume of the well-known Die Differential- und Integralgle- 
ichungen der Mechanik und Physik, edited by P. Frank with the assistance 
of H. Faxen, R. Fiirth, Th. von Karman, R. von Mises, Fr. Noether, C. 
W. Oseen, A. Sommerfeld, and E. Trefftz, Braunschweig, Vieweg, 1927-35, 
whose second volume is a revised and enlarged edition of the second part 
of the celebrated 1910 Riemann-Weber treatise Die Differential- und Inte- 
gralgleichungen der Mechanik und Physik. Hence, Sommerfeld transforms 
(*) in a complex integral, and, to this end, on the basis of the relation 

23 Here, ( n (x) = ttx/2H((((_ 1 i 2 (x) (where Ha\x) is the Hankel function of the second 

kind), is the so-called Riccati-Bessel function, while £ n (x) = Cn{x) +x^(a;). 
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P n (cos6 ) = (—l) n P n (—cos9), which is valid for integral (and only for inte¬ 
gral) values of n, he first rewrites the series in (*) in the following form 


(*i) 


u = V(2n+l)(-l)"P n (-cosl9)^ 

n=o 


hence, he replaces n by a complex variable v and traces a loop A, around 
the real axis of the //-plane, in such a manner it surrounds all the points 
v = 0,1, 2, 3, n ,... in a clockwise direction. Over this loop, Sommerfeld 
takes the integral 


(*2) 


u = 


2v + 1 
2i sin vi r 


P u { — cos 6) 


Cn(kr) 

Uka) U 


which is obtained from the general term in (nq) by interchanging n and v, 
suppressing the factor (—l) n , and appending the denominator sini m. As 
in (Sommerfeld 1949, Chapter V, Appendix II), P v does not stand for the 
Legendre polynomial, but for the hypergeometric function 

P u (x) = F — v,v + 1,1, —) 

which is identical with the Legendre polynomial only for integral v. In do¬ 
ing so, Sommerfeld resembles 1888 Pincherle’s method. Then, Sommerfeld 
suitably deforms, in many ways, the path A taking into account poles and 
zeros of the integrand of (* 2 ) and applying Cauchy’s residue theorem, hence 
discussing the related consequences in relation to (* 2 ) until np to obtain the 
following expression 


(*3) 





2v + l 
2i sin vn 


P u ( — cos 9) 


C n(kr) 

Vn(ka) 


where, in a neighborhood of the m-th pole of the integrand of (* 2 ), r h — 
{d^v/dv) v=Vm and $, u (ka) — (v — v m )rj„(ka), if one considers the zeros of 
$,„(ka) = 0 for v = n 0 , v \,... as the poles of the integrand of (* 2 ). Now, 
except for a sign and for a constant factor, the integral (* 2 ) is identical with 
the solution (=iq) of the sphere problem. Hence the series (*3) also represents 
the solution of the sphere problem, and suppressing the immaterial constant 
factor, Sommerfeld finally writes 


(*4) 


+ 1 D , Q ^ (u(kr) 

U = — - —Pn{— COS 9)- 


sm vi r 


T] u {ka) ’ 
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so that the passage from the series (*i), which is summed over integral n, to 
the series (*4), which is summed over the complex u, is obtained by forming 
residues in a complex integral twice. This, finally, explains why De Alfaro 
and Regge, in (De Alfaro & Regge 1965, Chapter 9, Section 9.3) (see also 
(De Alfaro et al. 1973, Chapter 1, Section 7)) speak of Watson-Sommerfeld 
transform (in short, WS-transform or SW-transform) , even if, in this paper, 
we have tried to deepen what have been the real historical roots of this 
transform, identifying these in the previous works of E.W. Barnes and S. 
Pincherle. 

In conclusion, in this last section we have only pointed out some lit¬ 
tle known historical perspectives which underlie the introduction of WS- 
transform, identifying the early roots upon which it relies. As we have seen in 
the previous Section 1, after the pioneering use of the WS-transform involving 
complex angular momenta for studying resonance phenomena of elementary 
particles mainly motivated to prove the validity of the Mandelstam represen¬ 
tation for the potential scattering of two spinless particles for a certain class 
of generalized Yukawa potentials, it was possible to work out the notions 
of Regge pole and Regge trajectory, from which Gabriele Veneziano started 
his as many pioneering work of 1968 (see (Veneziano 1968)) in which reso¬ 
nance dual models are introduced, for the first time, in fundamental physics. 
Veneziano’s model basically begins with the proposing of a quite simple ex¬ 
pression for the relativistic scattering amplitude (today known as Veneziano 
amplitude ) that obeys the requirements of Regge asymptotics and crossing 
symmetry in the case of linearly-rising trajectories, containing automatically 
Regge poles in families of parallel trajectories with residue in definite ratios, 
which furthermore satisfies the conditions of super convergence (S. Fubini) 
and exhibits, in a nice fashion, the duality between Regge poles and reso¬ 
nances in the scattering amplitude, essentially based on symmetry properties 
of Euler beta function (see (Green et al. 1988, Chapter 1, Sections 1.1 and 
1.2)). This is the incipit of a coherent S matrix theory of strong interactions 
which then will led to resonance dual models in high-energy particle physics, 
so, as is well known, marking the rising of string theory. Indeed, following 
(Di Vecchia 2008), once that Regge discovered that the Schrodinger equation 
allowed continuation in angular momentum for complex values, and linked 
resonances with different spin, one of the basic ideas that led to the construc¬ 
tion of an S matrix was that it should include resonances at low energy and at 
the same time give Regge behavior at high energy. But the two contributions 
of the resonances and of the Regge poles should not be added because this 
would imply double counting. This was called Dolen-Horn-Schmidt duality, 
discovered around 1967, while another idea (of about 1968) that helped in 
the construction of an S matrix was the so-called planar duality that was 
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visualized by associating to a certain process a duality diagram where each 
meson was described by two lines representing the quark and the antiquark. 
Finally, also the requirement of crossing symmetry played a very important 
role, and starting from these ideas Veneziano was finally able to construct 
an S matrix for the scattering of four mesons that, at the same time, had an 
infinite number of zero width resonances lying on linearly-rising Regge tra¬ 
jectories and Regge behaviour at high energy. In any case, after the Regge 
1959 seminal paper, the WS-transform was extensively used in field theory 
as, for instance, done in (Schwarz 1973, Section 1), where, on the basis of 
an axiomatic scheme (amongst which the assumption that all mesons and 
baryons, together with the related poles, lie on a Regge trajectory), it is 
considered the following dispersion relation 


A(s,t) 


X 


71=0 


Rnjt ) 

n — a(s) 




to which the scattering amplitude must satisfy, together with its WS-transform, 
namely 


A(s,t) = — / m 


r(-n) 

n — a(s) 


whose WS-contour C will be suitably deformed to avoid the pole at n = a(s). 
Through a suitable choice of the function R(t), Schwarz derives Veneziano 
amplitude just by means of such a method based on WS-transform. Likewise 
to the paper (Mandelstam 1974), where, although not explicitly mentioned, 
the WS-transform plays a fundamental role since the beginnings. Further¬ 
more, following (Cushing 1990, Chapter 5, Section 5.4), the asymptotic be¬ 
havior of f(k, 6) as cos 9 —* oo, was settled for potential scattering by Regge 
in 1959. Although a beautiful piece of classical mathematical analysis in 
its own right, this paper was to have its greatest impact on the S matrix 
program for the conjectures to which it would lead. Only (Collins 1977, 
Chapter 2, Sections 2.7, 2.9 and 2.10; Chapter 4, Section 4.6; Chapter 7, 
Section 7.1; Chapter 9, Section 9.3) has devoted the right attention to the 
WS-transform, there called Mandelstam-Watson-Sommerfeld transform, af¬ 
ter the improvement of the original WS-transform by Mandelstam himself in 

1962, recalling as well the links between it and the Mellin transform; in this 
regard, see also (Collins & Squire 1968, Chapter II, Section II.9) and (Omnes 
& Froissart 1963, Chapter 4, Sections 4-2 and 4-3), where interesting histori¬ 
cal notes on Sonnnerfeld original work have also been inserted. In (Frautschi 

1963, Chapter X), as regard Watson-Sonnnerfeld transformation, it is then 
said that this method had a history stretching over several decades and had 
been used to study rainbows, propagation of radio waves around the Earth, 
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and scattering from various potentials, the author specifying that the related 
transformation was used by Poincare and Nicholson in 1910 in connection 
with the bending of electromagnetic waves by a sphere, the transformation 
having been introduced in its present form by Watson in 1918 and later res¬ 
urrected by Sonnnerfeld. Regge’s original contribution lays in understanding 
the special features of complex angular momenta for scattering from super¬ 
positions of Yukawa potentials - the type of potential believed relevant to 
relativistic scattering - and calling these features to the attention of high- 
energy physicists. In (Frautschi 1963, Chapters X and XIII), besides further 
and interesting its applications in high-energy physics, it is also quoted who 
has later improved and extended WS-transformation. 

Finally, in this paper, we have pointed out a single historical aspect un¬ 
derling one the most important works of 20th century physics, that is to say, 
the 1959 seminal paper of Tullio Regge. Following (Schwarz 1975), Regge’s 
pioneering work of 1959 provided, amongst other things, a remarkable clas¬ 
sification of hadrons by means of the so-called Regge’s trajectories. Further¬ 
more, as has been seen above, this Regge’s work has mainly a technically 
fashion in which new formal methods and techniques have been introduced, 
besides to have introduced new pioneering and fruitful ideas and notions 
having physical nature. Among these techniques, just the so-called Watson- 
Sommerfeld transform, on which we have focussed our historical sight with 
the main purpose to identify the early origins and sources upon which rely 
such a complex transform. In pursuing this aim, we have found early sources 
of this formal method in some works dating back the early 1900s and mainly 
due to Poincare, Nicholson, Macdonald, and Barnes, until to reach to descry 
some prolegomena, or protohistoric origins, of the WS-method in some previ¬ 
ous works of Pincherle of the late 1880s. Because of the notable importance 
played by this transform in the formal framework of fundamental physics, we 
think that this our preliminary historical work, albeit centered on a single 
methodological aspect, is not altogether superfluous if nothing else for the 
simple fact that such a historical clarifications have not been treated in a 
whole and deep manner if not fully neglected. Above all, we have found a 
very few historical accounts on WS-transform, which we have seen to be a 
crucial technique laying out at the early foundational bases of resonance dual 
models of elementary particles. 
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